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ABSTRACT 

A  computer  model  which  simulates  Naval  Gunfire  Support 
of  Small  Independent  Action  Forces  (SIAF)  in  an  insurgent- 
contested  area  is  presented.  The  mechanics  of  the  model  are 
described  in  sufficient  detail  to  allow  understanding  by  any 
user  with  a  basic  understanding  of  Naval  Gunfire  Support  and 
computer  simulation  techniques.  A  detailed  guide  to  pre¬ 
paring  the  input  data  deck  required  for  simulation  of  a  user 
defined  problem  is  included. 

The  problem  definition  and  results  of  a  partial  fea¬ 
sibility  study  of  the  Seabased  Tactical  Deterrent  Force 
(STDF)  concept  using  the  simulation  model  are  presented  as 
a  means  of  providing  some  insight  into  the  feasibility  of 
the  STDF  concept  and  as  an  example  of  the  use  of  the  simula¬ 
tion  model. 
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I.  INTRODUCTION 


In  these  days  of  reduced  U.S.  Military  manpower  and 
hardware,  a  great  deal  of  attention  is  being  focused  on 
developing  the  optimal  military  posture  for  future  opera¬ 
tions.  One  of  the  main  factors  being  considered  in  devel¬ 
oping  this  posture  is  the  United  States*  announced  desire 
to  reduce  the  number  of  U.S.  troops  and  troop  support  fa¬ 
cilities  on  foreign  soil  while  still  maintaining  the  ability 
to  respond  with  an  adequate  amount  of  military  power  to  meet 
commitments  to  allied  countries.  Thus,  a  foreseeable  pos¬ 
sibility  is  that  the  U.S.  might  be  called  upon  for  emergency 
military-  aid  to  a  country  in  which  it  has  no  land-based 
troops,  aircraft,  or  artillery,  as  was  the  case  in  the 
Lebanon  and  Dominican  Republic  crises  of  the  recent  past. 

One  of  the  concepts  of  future  military  posture  which  is 
under  consideration  is  that  of  the-  Seabased  Tactical  Deter¬ 
rent  Force  (STDF) .  Integral  to  the  STDF  concept  is  the  con¬ 
cept  of  the  Small  Independent  Action  Force  (SIAF) .  Each  of 
these  concepts  is  discussed  below. 

A.  SEABASED  TACTICAL  DETERRENT  FORCE  ' 

♦  ’ 

The  STDF  would  be  required  .to  have  a  multitude  of  cap¬ 
abilities.  Without  defining  all  of  the  capabilities  re¬ 
quired,  it  can  be  summarily  stated  that  the  STDF  would  have 
to  combine  the  existing  capabilities  of  an  Attack  Carrier 
Striking  Force,  an  Amphibious  Task  Force,  a  Naval  Gunfire 
Support  Force  and  an  Underway  Replenishment  and  Logistics 
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Force,  while  providing  for  its  own  defense  as  well.  One  of 
the  most  important  requirements  of  the  STDF,  and  the  one 
which  emphatically  distinguishes  it  from  other  alternative 
concepts,  would  be  the  ability  to  strategically  interpose 
combat  forces  in  most  areas  in  which  conflicts  might  arise 
which  would  seriously  jeopardize  vital  U.S.  interests.  This 
is  not  a  new  requirement,  but  it  would  enjoy  renewed  and 
expanded  importance  under  the  STDF  concept.  The  convention¬ 
al  approach  to  exercising  this  capability  which  has  been 
used  many  times  in  the  past  has  been  to  land  large  numbers 
of  troops  ashore  by  an  amphibious  landing  to  conduct  massive 
operations  to  achieve  a  military  objective.  But  with 
changing  U.S.  foreign  policy,  such  massive  operations  ap¬ 
pear  less  likely.  Also,  our  experience  in  Southeast  Asia 
has  inspired  new  technology  and  concepts  of  amphibious  war¬ 
fare.  In  the  face  of  an  insurgent  force  employing  guerrilla 
warfare  or  similar  tactics,  massive  amphibious  operations 
are  rarely  practical  or  productive.  The  effectiveness  of 
small  reconnaissance  and  combat  patrols,  such  as  the  Small 
Independent  Action  Force  (SIAF)  patrols,  in  seeking  out  the 
insurgent  forces  has  often  been  demonstrated.  Thus,  the 
insertion  and  support  of  SIAF  patrols  becomes  a  capability 
which  the  STDF  must  have. 

B.  SMALL  INDEPENDENT  ACTION  FORCE  (SIAF)  PATROLS 

At  this  point  in  the  discussion  it  is  presumed  that  it 
has  been  determined  that  the  STDF  concept  is  a  viable  means 
of  interposing  U.S.  military  power,  and  the  ability  to  insert 
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and  support  SIAF  patrols  is  organic  to  the  STDF.  The  fol¬ 
lowing  question  then  arises:  ’’What  type  of  operations  by 
the  SIAF  patrols  would  be  most  effective  in  a  limited  war 
situation  if  one  arises?”  Some  high-ranking  military  of¬ 
ficers  have  expressed  the  opinion  that  it  might  be  tactic¬ 
ally  advantageous  in  many  areas  to  saturate  an  area  of 
limited  size  with  SIAF  patrols.  These  patrols  would  be 
dispersed  throughout  the  area  with  the  objective  to  either: 

1.  Reconnoiter  an  area  where  enemy  strength  and  degree 
of  control  is  virtually  unknown,  or 

2.  Seek  out  known  or  suspected  enemy  forces. 

In  either  case,  each  individual  patrol  should  typically 
be  of  short  duration  (a  few  days) ,  small  in  numbers  (4  to  18 
men) ,  dependent  only  on  the  STDF  for  logistics  and  external 
fire  support,  self-sustaining  while  actually  on  patrol,  and 
capable  of  rapid  and  reliable  communication  with  the  STDF. 

The  primary  mission  envisioned  for  the  SIAF  patrols  would 
be  one  which  is  essentially  passive.  The  patrols  would 
seek  out  enemy  personnel,  base  camps,  support  facilities, 
etc. ,  but  rather  than  actually  engaging  in  combat  with  such 
targets,  the  patrols  would  call  in  air  strikes  or  naval  gun¬ 
fire  from  the  STDF.  Thus,  the  mission  of  the  patrols  is 
primarily  reconnaissance/ forward  observer. 

C.  OBJECTIVES  OF  THIS  THESIS 

In  this  paper,  a  small  segment  of  the  overall  operations 
analysis  effort  required  to  investigate  the  STDF  and  SIAF 
concepts  is  isolated  and  cu..ducted.  Specifically,  the  subject 
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under  study  is  Naval  Gunfire  Support  (NGFS) ,  originating 
from  the  STDF,  in  support  of  SIAF’s  ashore  in  an  insurgent- 
contested  area.  The  objective  of  the  paper  is  twofold: 

1.  Development  of  Computer  Simulation  Model 

A  computer  simulation  model  has  been  developed 
which  simulates  the  actions  of  the  STDF  in  providing  NGFS 
to  SIAF’s  ashore.  The  model  is  designed  primarily  for  in¬ 
vestigating  the  queueing  aspects  of  providing  NGFS.  The 
main  focus  of  attention  is  the  time  interval  between  the 
time  a  SIAF  calls  for  a  mission  to  be  fired  and  the  time 
a  ship  commences  firing  the  mission.  No  ascertainment  is 
made  of  the  actual  results  achieved  when  the  mission  is 
fired..  The  model  is  very  flexible  in  design ,  allowing  the 
simulation  of  a  wide  variety  of  situations  defined  by  user- 
specified  inputs. 

2.  Partial  STDF  Feasibility  Study 

After  the  simulation  model  had  been  developed,  a 
typical  situation  in  an  insurgent -contested  area  was  hy¬ 
pothesized,  and  the  model  was  exercised  to  determine  the 

optimum  number  and  mix  of  various  type  *  of  NGFS  ships  which 

• 

would  be  required  to  adequately  support  the  SIAF’s  ashore. 
Since  situations  in  different  area^  certainly  vary  greatly 
and  may  be  more  or  less  severe  than  the  hypothesized  situ¬ 
ation,  and  since  the  simulation  is  restricted  solely  to  NGFS 
and  no  other  organic  capabilities  of  the  STDF  (such  as  air 
strikes) ,  the  results  of  the  exercise  cannot  be  considered 
to  provide  any  definite  conclusion  concerning  the  overall 
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feasibility  of  the  STDF  concept.  However,  the  exercise 
provides  an  example  of  the  use  of  the  simulation  model  for 
other  users  to  follow  when  a  specific  real-life  situation 
is  under  investigation,  and  it  also  provides  some  insight 
into  the  feasibility  of  a  STDF  providing  NGFS  in  support 
of  SIAF's  ashore  in  an  insurgent- contested  area. 

D.  ORGANIZATION  OF  THIS  THESIS 

The  development  and  use  of  the  model  are  discussed  in 
this  thesis.  In  the  next  section,  the  systems  to  be  simu¬ 
lated  are  described  to  the  degree  of  detail  necessary  to 
construct  the  simulation  model.  Then,  the  mechanics  of  the 
model  are  discussed  in  sufficient  detail  to  allow  any  user 
with  a  basic  knowledge  of  NGFS  and  simulation  techniques 
to  understand  how  the  computer  program  works.  Next,  the 
method  and  results  of  a  partial  STDF  feasibility  study  are 
presented  and  conclusions  based  on  the  results  are  dis¬ 
cussed.  Finally,  refinements  and  extensions  of  the  model 
are  discussed. 
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II.  tHE  system 


The  primary  system  which  is  simulated  can  be  identified 
as  Naval  Gunfire  Support.  In  order  to  construct  a  realistic 
and  useful  simulation  model,  at  least  two  other  systems  must 
be  included  in  the  simulation  model.  These  are  refueling 
and  rearming  of  ships  at  sea  and  the  conduct  of  SIAF  patrols 
ashore.  Each  of  these  systems  is  described  in  this  section 
of  the  paper. 

A.  NAVAL  GUNFIRE  SUPPORT 
1.  NGFS  Ship  Types 

Any  U.S.  Naval  ship  with  installed  guns  is  considered 
capable  of  providing  some  degree  of  gunfire  support  to  troops 
ashore.  However,  from  among  the  many  types  of  U.S.  Naval 
ships  with  installed  guns,  the  only  ones. which  can  be  con¬ 
sidered  to  have  NGFS  as  one  of  their  primary  duties  are 
cruisers  (CA,  CL,  CLG) ,  destroyers  (DD,  DL,  DDG,  DLG) ,  and 
battleships  (BB) .  In  recent  years,  certain  types  of  ships, 
such  as  the  Inshore  Fire  Support  Ship  (IFS)  and  the  Medium 
Rocket  Launching  Ship  (LSMR) ,  have  been  built  or  converted 
with  the  specific  primary  mission  of  providing  NGFS.  How¬ 
ever,  these  ships  normally  do  not  have  adequate  cruising 
ranges  and  speeds  to  be  considered  as  permanent  elements  of 
a  STDF.  There  are  no  battleships  in  commission  at  the  . 
present  time.  Therefore,  the  only  types  of  NGFS  ships  con¬ 
sidered  in  the  model  are  cruisers  and  destroyers. 
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2.  Installed  Guns 


All  cruisers  and  destroyers  have  more  than  one 
type  of  gun  installed,  but  each  has  one  type  which  is  con¬ 
sidered  to  be  it’s  "main  battery"  for  NGFS.  For  example, 
a  heavy  cruiser  (CA)  normally  has  6  8"/55  guns,  10  5"/38 
guns,  and  8  3"/50  guns.  However,  the  8"/55  guns  are  con¬ 
sidered  to  be  the  main  battery  and  would  normally  be  the 
only  ones  used  for  NGFS,  while  the  other  guns  are  primarily 
used  as  a  backup  and  for  self-defense  in  surface-to-air 
and  surface-to-surface  engagements.  In  the  model,  each 
ship  is  considered  to  use  only  its  main  battery  for  NGFS. 
The  types  of  ships  considered  in  the  model  and  the  main 
battery  gun  type  of  each  [Reference  1]  are  listed  below: 


Cruisers 

Destroyers 

CA 

8"/55 

DD  (Fram  I/II) 

5"/38 

CL 

6"/47 • 

DD  (931  Class) 

5*754 

CLG 

6"/47 

DDG  ' 

5"/54 

DLG 

5*’/54 

DL 

5*754 

Throughout  this  paper,  a  CA  is  referred  to  as  an 
8"/55  ship,  a  CL  or  CLG  as  a  6"/47^s.hip,  etc. 

3.  Firing  Procedures 

Procedures  for  calling  for  and  providing  NGFS  vary 
with  the  specific  situation  and  local  instructions.  In  some 
situations  ships  might  simply  proceed  to  a  firing,  position 
and  fire  upon  a  predetermined  area  without  the  assistance  of 
spotters  ashore.  This  is  called  harrassment  and  interdiction 


18 


(H§I)  fire,  but  this  type  of  fire  is  not  envisioned  in  the 
STDF/SIAF  concept.  The  type  of  firing  which  is  most  likely 
to  occur  in  the  STDF/SIAF  situations  can  be  considered  in 
almost  all  cases  to  adhere  to  the  following  chronological 
sequence  of  events: 

a.  Call  for  Fire 

A  SIAF  patrol  detects  a  target  and  radios  a 
request  for  NGFS  to  the  STDF.  This  is  referred  to  as  a 
'’call  for  fire”,  and  will  be  abbreviated  as  CFF  throughout 
this  paper. 

b.  Ship  Assigned 

The  STDF  commander  or  his  delegated  subordinate 
determines  which  ship  from  among  those  available,  if  any  are 
available,  is  most  appropriate  to  fire  the  mission  and  as¬ 
signs  that  ship.  This  assignment  is  based  on  the  following 
factors : 

(1)  Range  to  the  target, 

(2)  type  of  the  target, 

(3)  Distance  of  each  ship  from  the  firing  posi¬ 
tion. 

c.  Assigned  Ship  Proceeds  to  Firing  Position  . 

The  firing  position  for  each  mission .is  considered 
to  be  the  point  along  the  "gunline"  which  is  closest  to  the 
target  (the  perpendicular  distance  from  the  gunline  to  the 
target) .  The  gunline  is  a  fictitious  line  in  the  ocean  run¬ 
ning  roughly  parallel  to  the  shoreline  at  a  specified  dis¬ 
tance  from  the  shore. 
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d.  Assigned  Ship  Sets  Up  for  Mission 

After  the  assigned  ship  has  reached  the  firing 
position,  it  fixes  its  navigational  position  and  readies 
its  gun  and  then  reports  to  the  SIAF  that  it  is  ready  to 
fire . 

e.  Fire  the  Mission 

The  assigned  ship  commences  firing  on  a  signal 
from  the  SIAF.  After  the  first  few  rounds  have  been  fired, 
the  SIAF  observes  the  fall  of  shot  and  radios  instructions 
to  the  ship  to  bring  fire  onto  the  target.  This  is  called 
"spotting".  When  fire  has  been  adjusted  onto  the  target, 
the  ship  continues  to  fire  at  a  rapid  rate  until  the  SIAF 
observes  that  the  target  has  been  destroyed  or  otherwise 
effectively  encountered.  The  ship  then  ceases  fire  and  re¬ 
turns  to  a  "ready"  status. 


Ships  assigned  to  NGFS  duties  but  not  actually  en¬ 
gaged  in  firing  would  normally  remain  in  the  vicinity  of  the 
gunline  in  order  to  be  able  to  respond  rapidly  when  assigned 
to  fire  a  mission.  The  specific  area  in  which  non- firing 
ships  would  operate  would  be  specified  by  the  STDF  commander. 
A  frequently  used  procedure  is  for  non- firing  ships  to  oper¬ 
ate  on  or  inside  of  a  race-track  type  waiting  pattern  at  a 
distan  t  from  the  gunline  determined  by  considering  the 
tradeoffs  between  rapid  response  time  and  the  need  to  keep 
non-firing  ships  out  of  range  of  suspected  enemy  shore  bat¬ 
teries  . 
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B.  REFUELING  AND  REARMING 

A  definite  limitation  on  the  availability  of  NGFS  ships 
is  imposed  by  fuel  and  ammunition  capacities  of  each  ship. 
Under  normal  steaming  conditions,  a  destroyer  with  a  full 
load  of  fuel  can  steam  for  3  to  4  days  before  it  will  be  re¬ 
quired  to  refuel,  while  a  cruiser  can  steam  5  to  6  days. 

This  assumes  that  refueling  is  required  when  the  amount  of 
fuel  on  board  reaches  501  of  ships  capacity,  a  fairly  com¬ 
mon  criterion.  To  determine  when  rearming  is  required,  the 
STDF  commander  would  normally  specify  a  minimum  on-board 
ammunition  level  to  which  ships  would  be  allowed  to  fall  be¬ 
fore  rearming.  This  level  would  be  a  function  of  the  type 
of  ship,  expected  rate  of  fire,  and  number  of  other  ships 
present. 

It  is  possible  that  refueling  and  rearming  requirements 
would  not  impose  a  limitation  on  NGFS  provided.  This  would 
be  the  case  when  other  NGFS  ships  ‘are  available  to  tempor¬ 
arily  replace  ships  which  are  required  to  refuel  or  rearm, 
or  when  the  duration  of  the  operation  is  short  enough  that 
rearming  and  refueling  should  not  be  necessary. 

C.  SIAF  PATROL  PROCEDURES  %  . 

A  very  detailed  computer  simulation  model  of  SIAF  patrol 
actions  has  been  developed  by  TRW,  Inc.  of  Los  Angeles  and 
is  described  in  detail  in  Reference  2.  This  model  can  be 
used  for  developing  optimal  patrol  procedures  and  for  deter¬ 
mining  expected  patrol  results  for  a  wide  variety  of  friendly 
and  enemy  situations.  In  the  model  described  in  this  thesis, 
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no  attempt  has  been  made  to  simulate  specific  patrol  actions 
with  the  exception  of  patrol  insertion,  calls  made  to  the 
STDF  for  NGFS,  and  patrol  extraction.  Under  the  STDF  con¬ 
cept,  a  SIAF  patrol  would  typically  adhere  to  the  following 
chronological  sequence  of  major  events: 

1.  Insertion 

The  patrol  is  inserted  into  its  assigned  subarea  of 
the  overall  objective  area.  This  insertion  would  normally 
be  accomplished  by  helicopter  from  an  LHA  assigned  to  the 
STDF,  but  could  conceivably  be  accomplished  by  other  means 
such  as  paradrop  or  amphibious  landing.  Order  and  time  of 
insertion  into  the  various  subareas  would  be  chosen  such  as 
to  yield  the  highest  probability  of  achieving  a  tactical  ad¬ 
vantage  by  surprising  and  confusing  the  enemy. 

2.  Patrol  Maneuvers 

The  SIAF  patrol  maneuvers  within  its  assigned  sub- 
area  attempting  to  locate  the  enemy  and  his  support  facil¬ 
ities  while  remaining  undetected  by  the  enemy.  When  targets 
are  located  the  patrol  calls  for  NGFS  or  air  strikes,  if 
available,  rather  than  actually  engaging  in  combat  with  the 
enemy.  The  patrol  would  continue  to  maneuver  within  its  as¬ 
signed  subarea  until  extracted  for  one  of  the  three  reasons 
indicated  below. 

3.  Extraction 

The  patrol  is  extracted  when  one  of  the  following 

occurs  * 

a.  A  predetermined  extraction  time  is  reached. 
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b.  The  patrol  leader  considers  that  optimal  area 
coverage  has  been  achieved. 

c.  Extraction  is  required  because  the  patrol's 
effectiveness  has  been  significantly  reduced  by  personnel 
casualties . 

Extraction  would  normally  be  accomplished  by  heli¬ 
copter. 

D.  THE  ENEMY  SITUATION 

Obviously,  the  enemy  situation  is  a  factor  which  plays 
a  large  role  in  any  study  of  NGFS  and  SIAF  patrols.  No 
study  could  ever  begin  without  a  reasonably  accurate  estima¬ 
tion  of  the  size  of  the  enemy  force,  its  locations,  its  cap¬ 
abilities,  and  its  expected  actions.  These  enemy  factors 
and  others  are  considered  in  great  detail  in  the  TRW  SIAF 
model.  They  are  considered  in  the  model  described  in  this 
paper  in  the  following  ways: 

1 .  Subarea  Division 

Division  of  the  objective  area  into  patrol  sub- 
areas  is  based  on  terrain  features,  expected  enemy  locations, 
and  expected  number  and  capabilities  of  enemy  in  each  sub- 
area. 

2 .  Call-for-Fire  Rate 

The  real  cornerstone  on  which  the  model  depends  for 
accurate  estimates  of  the  queueing  aspects  of  NGFS  is  the 
expected  rate  of  calls  for  fire  in  each  subarea.  This  rate 
is,  of  course,  primarily  dependent  on  the  enemy  density  in 
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each  subarea.  For  a  specific  situation,  the  CFF  rate  for 
a  particular  subarea  can  be  closely  estimated  by  the  TRW 
SIAF  model  or  can  be  predicted  based  on  available  intelli¬ 
gence  . 

E.  ASSUMPTIONS 

There  are  several  general  assumptions  which  are  neces¬ 
sary  in  modeling  NGFS  within  the  scope  of  this  paper.  These 
are  discussed  below.  While  each  of  these  detracts  from  the 
realism  of  the  model,  it  is  considered  that  none  of  them 
significantly  detracts  from  the  accuracy  of  queueing  re¬ 
sults.  Other  more  specific  assumptions  will  be  discussed 
as  they  occur  in  the  description  of  the  model  in  Section 
III. 

1.  SIAF  patrols  remain  within  their  assigned  subarea 
throughout  the  time  interval  between  insertion  and  extrac¬ 
tion,  and  only  detect  targets  in  their  assigned  subarea. 

2.  The  enemy  situation  at  the  start  of  the  operation 
in  each  subarea  remains  the  same  throughout  the  operation. 

3.  The  enemy  does  not  launch  attacks  against  NGFS 
ships,  or  if  he  does  the  attacks  are  not  successful  enough 
to  significantly  reduce  the  operational  capability  of  the 
ship. 

4.  SIAF  patrols  do  not  come  under  enemy  attack,  or  if 
they  do  the  attack  is  such  that  emergency  NGFS  is  not  re¬ 
quired  and  the  patrol's  capability  is  not  seriously  reduced. 

5.  A  SIAF  patrol  remains  stationary  at  the  same  point 
from  the  time  it  has  called  for  a  mission  until  the  mission 
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is  complete.  This  leads  to  the  further  assumption  that 
# 

patrols  spot  only  one  target  at  a  time,  i.e.,  while  one 
mission  is  in  progress  the  patrol  does  not  detect  any 
other  targets. 

6.  When  a  SIAF  patrol  is  extracted  from  its  assigned 
subarea,  another  patrol  is  not  inserted  into  the  area  to 
replace  the  extracted  patrol. 

7.  SIAF  patrols  do  not  sleep;  at  least,  no  element  of 
the  model  specifically  simulates  a  period  of  time  during 
which  patrols  are  sleeping.  However,  it  can  be  considered 
that  unusually  long  interarrival  times  between  missions  in 
a  particular  subarea  could  be  caused  in  some  instances  by 

a  period  of  time  during  which  the  patrol  in  that  subarea  is 
sleeping. 


III.  THE  MODEL 


This  section  describes  the  model  which  has  been  devel¬ 
oped  to  simulate  the  systems  described  in  Section  II.  The 
simulation  program  is  written  in  FORTRAN  IV.  It  is  not  the 
intent  of  this  section  to  describe  the  entire  program  step- 
by-step;  but  rather  to  describe  the  mechanics  of  the  program 
in  sufficient  detail  to  allow  any  user  with  a  basic  knowledge 
of  NGFS  and  computer  simulation  techniques  to  understand  what 
is  being  simulated  and  to  be  able  to  use  the  model  according 
to  his  own  input  parameters ,  or  to  adapt  certain  sections  of 
the  program,  if  desired.  A  generalized  flow  chart  of  the 
program  is  contained  in  Appendix  A  and  the  complete  program 
listing  is  contained  in  Appendix  3.  Numerous  comment  cards 
have  been  included  in  the  program  listing  to  facilitate  un¬ 
derstanding. 

A.  PRELIMINARY  DISCUSSION 

1.  Variable  Names 

Mnemonic  variable  names  are  used  throughout  the  pro¬ 
gram  to  facilitate  association  between  the  names  and  the 
quantities  which  they  represent. 

2.  Measures 

All  times  in  the  program  are  measured  in  whole 
minutes  and  distances  are  measured  to  the  nearest  whole  yard, 
unless  otherwise  indicated.  In  view  of  these  measures,  all 
variables  beginning  with  A  through  W  are  declared  to  be  4- 
byte  integers.  Whenever  a  calculation  of  an  integer-valued 


time  quantity  is  performed,  0.5  is  aided  to  the  result,  re¬ 


sulting  in  a  truncated  integer  accurate  within  ±.5  minutes. 


3.  Input  Data-General 

Almost  all  numerical  quantities  used  in  program  cal¬ 
culations  are  read  in  at  the  beginning  of  the  program  from  a 
data  deck  prepared  by  the  user.  The  definition  and  use  of 
each  input  variable  is  discussed  as  it  appears  in  the  model 
description  which  follows.  A  detailed  guide  to  assembling 
the  data  deck  is  included  in  Section  III.J. 

4.  Random  Numbers 

There  are  14  points  in  the  program  at  which  random 
numbers  are  used  for  Monte  Carlo -type  comparisons  or  calcula¬ 
tions.  Random  number  streams  at  each  point  are  independent 
of  all  other  points  in  the  program.  The  random  number  gen¬ 
erator  is  initialized  by  specifying  any  14  odd  five-digit 
numbers  as  input  variables  IR1  through  IR14.  KR  is  set  equal 
to  65539  at  the  beginning  of  the  program.  Random  numbers 
are  then  produced  by  the  following  calculation  at  each  point 
a  random  number  is  required  (IR1  is  used  for  example) : 


IR1  *  IR1  *  KR 

XRN  »  0.5  +  FLOAT (IR1)  *  2.328306E-10 

XRN  will  be  a  floating  point  random  number  uniformly 
distributed  between  0.0  and  1.0  [Reference  3]. 

5.  Ship  Identification 

Throughout  the  program  ships  are  referred  to  by 
their  type  (T)  and  ship  number  (S) .  The  four  types  of  ship 
and  corresponding  values  of  T  are: 
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1  ■  8*755  ship  (cruiser) 

2  -  6*747  ship  (cruiser) 

3  *  5*'/54  ship  (destroyer) 

4  ■  5*738  ship  (destroyer) . 

A  maximum  number  of  9  ships  of  each  type  are  pro¬ 
vided  for,  or  a  maximum  of  36  ships.  The  number  of  each 
type  of  ship  ''present  in  the  problem"  is  specified  by  the 
input  variables  NR8S5,  NR647,  NR554,  and  NR538.  "Present 
in  the  problem"  refers  to  the  number  of  ships  present  at 
the  start  of  the  problem.  The  names  of  quantities  which 
are  related  to  a  particular  ship  are  subscripted  by  the 
ship  number  and  type  in  that  order,  i.e.,  UTIL  (5,2)  is  the 
cumulative  total  minutes  that  6"/47  ship  #5  has  been  uti¬ 
lized  for  NGFS. 

6.  Uniform  Distributions 

There  are  several  points  in  the  program  where  ob¬ 
servations  of  a  random  variable  from  a  uniform  distribution 
are  made.  In  each  case,  the  observations  are  generated  by 
use  of  the  following  formula: 

RV  =  MIN  +  XRN  *  (MAX  -  MIN)  . 

where  RV  is  the  random  variable  being  observed,  XRN  is  a  ran 
dom  number  between  0.0  and  1.0,  and  MIN  and  MAX  are  the  pa¬ 
rameters  of  the  uniform  distribution. 

7.  Number  of  Iterations 

The  user  specifies  the  number  of  iterations  of  the 
problem  which  he  desires  to  use  to  obtain  averaged  results. 
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This  number  is  specified  as  the  input  variable  NRIT.  In 
several  trial  runs  of  the  program,  NRIT  =  10  was  shown  to 
be  sufficient  for  obtaining  reliable  results  with  a  rela¬ 
tively  narrow  951  confidence  interval.  However,  increasing 
NRIT  will  reduce  the  width  of  the  confidence  interval  even 
further. 

8.  Specification  of  Card  Reader  and  Printer  Logical 

Imit  Numbers 

All  READ  and  WRITE  statements  in  the  program  use  the 
values  of  the  variables  R  and  W  as  the  card  reader  and  prin¬ 
ter  logical  unit  numbers,  respectively.  As  written,  the 
program  specifies  R  =  5  and  W  =  6  at  the  beginning  by  the 
statement  DATA  R,W/5,6/.  If  the  card  reader  and  printer 
logical  unit  numbers  at  a  computer  facility  on  which  the 
program  is  to  be  run  are  not  5  and  6,  respectively,  then  the 
above  DATA  statement  must  be  replaced  by  one  which  specifies 
th?  appropriate  logical  unit  numbers. 

B.  PROBLEM  AREA 

the  first  step  in  developing  data  for  the  model  is  to 
define  the  problem  area  in  terms  of  a  coordinate  system 
which  can  be  specified  to  the  computer. 

1.  Objective  Area  COA) 

The  OA  is  the  overall  area  ashore  in  which  SIAF 
patrols  are  conducted.  The  OA  is  divided  into  NRAREA  patrol 
subareas,  where  NRAREA  is  an  input  variable  specifying  the 
number  of  subareas.  NRAREA  must  be  between  2  and  50.  These 
subareas  must  be  rectangular  and  the  boundaries  of  each  sub- 
area  must  be  parallel  or  perpendicular  to  the  boundaries  of 
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all  other  subareas,  and  the  northern  and  southern  boundaries 
of  each  subarea  must  be  parallel  to  the  gunline.  For  pur¬ 
poses  of  uniformity  it  is  considered  that,  all  subarea  bound¬ 
aries  run  north  to  south  and  east  to  west.  The  western  edge 
of  the  westernmost  subarea  is  considered  to  be  the  y-axis  of 
an  x-y  coordinate  system  (see  Section  III.B.3  for  exception). 
The  distance  of  the  eastern  ana  western  edges  of  each  sub- 
area  from  the  y-axis  are  the  input  variables  EAST(NR)  and 
WEST(NR),  where  NR  is  the  subarea  number.  Subareas  may  be 
numbered  in  any  way  desired  by  the  user. 

2.  Gunline 

The  gunline  is  a  fictitious  straight  line  in  the 
ocean,  parallel  to  a  straight  line  oriented  east-west  which 
approximates  the  shore  line.  See  Section  II.A.3.C  for  a 
definition  of  the  gunline.  The  gunline  is  considered  to  be 
the  x-axis  of  the  x-y  coordinate  system.  The  distance  of 
the  northern  and  southern  boundaries  of  each  subarea  from 
the  gunline  are  the  y-coordinates  of  the  subarea,  specified 
by  the  input  variables  NORTH (NR)  and  SOUTH (NR).  The  y-co- 
ordinate  of  the  northern  edge  of  the  northernmost  subarea 
must  be  greater  than  the  maximum  range  of  the  longest  gun 
present .  * 

3.  Waiting  Pattern 

The  area  in  which  ships  operate  while  not  assigned 
to  a  mission  is  defined  by  the  us«t  by  appropriate  specifica¬ 
tion  of  values  of  the  following  input  variables: 

WTDIS  -  Perpendicular  distance  of  the  northern  edge 
of  the  waiting  pattern  from  the  gunline 
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NSWAIT  -  The  north-south  length  across  the  waiting 

pattern 

WTNW  and  WTNE  -  The  perpendicular  distance  of  the 
western  and  eastern  edges,  respectively,  of  the  waiting  pat¬ 
tern  from  the  y-axis.  (NOTE:  If  the  western  edge  of  the 
waiting  pattern  is  west  of  the  western  edge  of  the  OA,  then 
the  western  edge  of  the  waiting  pattern  should  be  used  as 
the  y-axis  (WTNW  *  0)). 

4.  Justification  of  Problem  Area  Definition 

It  is  realized  that  the  scheme  of  problem  area 
definition  described  above  does  not  exactly  conform  to  most 
real-life  situations.  For  example,  area  boundaries  would 
often  conform  to  rivers  or  treelines  and  not  necessarily  be 
straight  lines.  However,  straight  line  approximations  such 
as  used  in  this  model  should  not  significantly  effect  the 
queueing  results  if  input  variables  are  properly  specified. 

Figure  1  illustrates  a  typical  problem  area  defined 
in  the  coordinate  system  described  in  this  section. 

C.  EVENT  CHAIN  AND  MISSION  SCHEDULE 

The  two  principle  mechanisms  by  which  data  is  kept  track 
of  and  program  control  is  exercised  are  the  ’’Event  Chain”  and 
the  ’’Mission  Schedule". 

1.  Event  Chain 

The  simulation  uses  the  "next-event"  method  of  time 
advance.  At  the  beginning  of  the  problem  the  statements 
K  »  1  and  CLOCK  *  EVENT (K)  sets  the  problem  timer  to  the  time 
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Example  of  Area  Coordinates: 

NORTH (1)  =  20000 
SOUTH (1)  =  13000 
WEST(l)  =  0 
EAST (1)  =  8000 


\  t 

True  Problem 
North  North 


Figure  1.  Typical  Problem  Area  Map. 
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of  the  first  event.  This  event  is  then  directed  through  its 
appropriate  path  of  logic.  When  it  has  been  completely  pro¬ 
cessed  the  value  of  K  is  increased  by  1  and  the  time  is  ad¬ 
vanced  to  the  time  of  the  next  event.  This  process  is 
continued  until  an  event  with  a  time  of  9999999  is  the  next 
event  and  the  problem  is  terminated.  Such  an  event  is 
created  when  the  generated  time  of  the  next  call  for  fire 
in  a  subarea  is  greater  than  the  time  of  extraction  of  the 
patrol  from  that  subarea. 

a.  Event  Types 

There  are  four  types  of  events  considered  An 
the  program.  These  are  indicated  by  the  value  of  the  vari¬ 
able  EVTYP(Kj  as  follows: 

1:  A  call  for  fire  from  a  SIAF  patrol 

2:  A  NGFS  mission  completion 

3:  A  scheduled  refueling  of  a  NGFS  ship 

4:  A  completion  of  refueling,  rearming,  or 
combined  refueling/ rearming  of  an  NGFS  ship. 

b.  Event  Parameters 

The  Event  Chain  actually  consists  of  four  param¬ 
eters  associated  with  each  event.  These  are  stored  in  vector 
arrays  with  variable  names  as  follbws : 

EVENT (K) :  Time  of  K^1  event  (minutes  from  zero 

problem  time) 

EVTYP(K) :  Type  of  K*h  event  (values  IAW  Section 

IU.C.l.a) 

EVTMSN(K) :  If  EVTYP(K)  =  1  or  2,  then  EVTMSN(K) 
is  the  mission  number  of  the  event.  Mission  numbers  are 
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assigned  in  chronological  order  (see  next  subsection):  If 
EVTYP(K)  =  3  or  4,  then  EVTMSN(K)  =  0. 

REPSHP(K) :  If  EVTYP(K)  =  1  or  2,  then  REPSHP(K) 
-  0.  If  EVTYP(K)  =  3  or  4,  then  REPSHP(K)  is  a  two  digit 
number  representing  the  number  and  type  of  the  refueling 
and/or  rearming  ship,  i.e.,  REPSHP(K)  =  23  means  5*754  ship 
#2  is  either  scheduled  to  refuel  or  has  completed  refueling 
and/or  rearming. 

c.  Subroutine  VENTCH 

Each  time  a  new  event  is  generated  in  the  pro¬ 
gram  subroutine  VENTCH  is  called.  The  arguments  of  VENTCH 
are  the  four  parameters  discussed  above.  S/R  VENTCH  inserts 
the  new  event  into  its  proper  chronological  order  on  the 
event  chr?in,  and  increases  the  event  number  of  all  suceeding 
events  by  1. 

2 .  Mission  Schedule 

The  mission  schedule  represents  a  chronological 
listing  of  all  NGFS  missions  which  have  been  generated  dur¬ 
ing  the  problem. 

a.  Mission  Generation 

A  call  for  fire  represents  a  mission  arrival. 

In  view  of  assumption  5  in  Section  II. E.,  mission  interarri¬ 
val  times  are  represented  by  the  time  interval  between  the 
time  the  last  mission  in  a  particular  subarea  was  completed 
and  the  time  of  the  next  mission  arrival  in  that  subarea. 
Mission  interarrival  times  in  each  subarea  are  assumed  to  be 
exponentially  distributed  with  a  rate  which  is  indirectly 


34 


specified  by  the  user  in  the  following  manner:  For  each  sub- 
area.  the  minimum  and  maximum  expected  patrol  duration  (num¬ 
ber  of  days  from  insertion  time  to  the  nearest  tenth)  are 
specified  by  the  user  as  the  input  variables  XSHORT(NR)  and 
XLONG(NR).  The  program  calculates  expected  patrol  duration 
as 

XDURAT(NR)  =  (XSHORT(NR)  +  XLONG (NR) ) / 2 . 

For  each  subarea,  the  user  also  specifies  the  expected  number 
of  calls  for  firing  during  a  patrol  (i.e.,  during  the  time 
XDURAT(NR)),  as  the  input  variable  XCFF(NR) .  Then  the  pro¬ 
gram  calculates  the  call-for-fire  rate  for  each  subarea  as 

XCALRT(NR)  =  XCFF (NR)/ (XDURAT(NR)  *  1440) 

where  the  factor  1440  converts  days  to  minutes.  The  arrival 
rate  then  has  units  of  calls/minute. 

Exponentially  distributed  interarrival  times 
are  generated  as  follows  (Reference  4]: 

The  cumulative  distribution  function  for  the 
exponential  distribution  is 

F(x)  =  1  -  e"Xx 

where  x  is  an  interarrival  time  and  X  is  the  arrival  rate. 
Since  F(x)  is  a  continuous  probability  between  0.0  and  1.0 
it  can  be  replaced  by  a  random  number  between  0.0  and  1.0, 
say  r.  Then 

r  -  1  -  e  *x  =$>x  «  - (l/X)ln(l-r) 
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or  equivalently 


x  *  -(1/A)ln(r)  . 

Thus,  -ALOG(XRN)/XCALRT(NR)  is  an  exponentially  distributed 
interarrival  time  of  a  call  for  fire  in  subarea  NR.  There 
are  two  points  in  the  program  at  which  calls  for  fire  are 
generated. 

(1)  First  Mission.  At  the  start  of  the  problem, 
the  time  of  the  first  mission  to  be  called  for  in  each  sub- 
area  is  calculated  as 

FIRST (NR)  «  INSERT (NR)  -  ALOG(XRN)/XCALRT(NR)  +0.5 

where  INSERT (NR)  is  the  time  of  patrol  insertion  into  subarea 
NR  in  minutes  from  zero  problem  time.  INSERT (NR)  is  an  in¬ 
put  variable  specified  by  the  user.  Zero  problem  time  is 
defined  as  the  number  of  minutes  after  midnight  on  the  first 
day  of  the  problem  that  the  first  partrol  is  landed  in  its 
assigned  subarea.  It  is  specified  by  the  input  variable 
HHOUR.  The  value  of  HHOUR  is  used  in  the  program  only  to 
convert  times  of  calls  for  fire  and  mission  completion  to 
date -time -groups  in  the  Mission  History  printout. 

These  first  missions  are  then  placed  on 

»  ■ 

the  mission  schedule  by  setting  MSNSKD(I)  «  FIRST(I)  and 
arranging  the  elements  of  MSNSKD  in  chronological  order  by 
a  "ripplesort"  technique.  The  first  missions  are  then  placed 
on  the  Event  Chain  with  parameters  EVENT(K)  *  MSNSKD(I), 
EVTYP(K)  =  1,  EVTMSN(K)  *  I,  and  REPSHP(K)  =  0. 
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(2)  Next  Missions.  Each  time  a  mission  is  com 
pleted,  the  time  of  the  next  call  for  fire  in  that  subarea 
is  calculated  as 

NEXT  *  CLOCK  -  ALOG (XRN) /XCALRT (AREA(E) )  +0.5 

where  CLOCK  is  the  present  value  of  the  problem  time  and 
AREA(E)  is  the  subarea  in  which  the  mission  was  just  com¬ 
pleted.  This  mission  is  then  placed  on  the  mission  schedule 
by  a  routine  which  is  similar  to  VENTCH  but  which  is  not  a 
separate  subroutine  since  it  is  only  used  at  this  one  point 
in  the  program.  Next  missions  are  also  placed  on  the  Event 
Chain  with  appropriate  parameters, 
b.  Mission  Parameters 

As  with  the  Event  Chain,  the  Mission  Schedule 
actually  consists  of  19  parameters  which  are  assigned  to 
each  mission.  When  a  mission  is  first  generated  it  is  only 
assigned  2  parameters,  MSNSKD(I)  and  AREA(I) .  The  other  17 
parameters  are  assigned  at  various  points  in  the  program 
after  CLOCK  time  has  reached  the  time  of  the  call  for  fire 
for  the  Itl1  mission.  A  discussion  of  each  of  the  vector  ar¬ 
rays  in  which  these  19  parameters  are  contained  and  the  way 
in  which  each  is  assigned  follows: 

(1)  MSNSKDm  :  Time  of  call  for  fire  for  the 
1*^  mission  (minutes  from  zerp  problem  time) 

(2)  AREA(I) :  Number  of  the  subarea  in  which 
the  mission  occurs 

(3)  RANGE(I) :  Range  from  the  firing  position 
to  the  target  of  the  I1*1  mission  (the  target’s  y-coordinate) 


(4)  EWTAR(I) :  East-west  coordinate  of  the  tar¬ 
get  of  the  I*k  mission  (the  target's  x-coordinate) 

Although  target  densities  vary  in  different 
subareas,  targets  are  considered  to  be  uniformly  distributed 
within  each  subarea.  Therefore  the  target  coordinates  are 
calculated  as  uniformly  distributed  random  variables  between 
SOUTH (A)  and  NORTH (A)  for  RANGE (I),  and  WEST (A)  and  EAST (A) 
for  EWTAR(I),  where  A  *  AREA(I) . 

(5)  TYPE (I) :  Type  of  target  of  the  I*h  mission 

There  are  four  target  types  considered  in  the 

problem.  These  are: 

1:  Many  personnel  (20  or  more)  or  personnel 
widely  dispersed. 

2:  Fisw  personnel  not  widely  dispersed. 

3:  Heavy  materia1  objects  (large  bunkers,  tanks, 
heavy  artillery,  etc.). 

4:  Light  material  objects  (small  bunkers,  jeeps, 
etc.). 

The  probability,  of  each  type  of  target  (i.e., 
the  estimated  percentage  of  all  targets  in  the  subarea  which 
are  of  each  type)  for  each  subarea  is  specified  by  the  user 
by  the  input  variables  XPROBl(NR),  XPR0B2(NR),  XPR0B3(NR)  , 
and  XPR0B4(NR).  The  sum  of  these  four  variables  must  equal 
1.0  for  each  subarea.  When  a  call  for  fire  occurs,  a  random 
number  is  generated  and  compared  with  the  cumulative  target 
type  probability  distribution  for  the  area  in  which  the  mis¬ 
sion  occurs  to  determine  the  target  type. 


38 


(6)  MISTYP(I) :  Mission  type  of  the  I*h  mission. 
This  type  refers  to  the  types  of  ships  which  may  be  assigned 
to  fire  the  mission,  as  will  be  discussed  in  Section  III.D. 

(7)  TYPSHP(I) :  Type  of  ship  assigned  to  fire 
the  1*^  mission. 

(8)  SHIPNR(I) :  Number  of  the  ship  assigned  to 
fire  the  Ith  mission. 

The  ways  in  which  parameters  (7)  and  (8)  are 
assigned  are  indicated  in  Sections  III.E.  and  III.F. 

(9)  TRAVEL(I) :  The  time  interval  between  the 
time  a  ship  is  assigned  to  fire  the  I**1  mission  and  the  time  * 
the  ship  arrives  at  the  firing  position. 

(10)  SETUP (I) :  The  time  interval  between  the 
time  the  assigned  ship  arrives  at  the  firing  position  and  the 
time  the  first  round  is  fired. 

(11)  FIRE(I) :  The  time  interval  between  the 
time  the  first  round  is  fired  and  the  time  the  last  round 
is  fired. 

(12)  SERVC(I)  :  Total  service  time  of  the  I**1 
mission.  SERVC(I)  =  TRAVEL (I )  +  SETUP(I)  +  FIRE(I). 

(13)  AMMEXP(I) :  Number  of  rounds  of  ammunition 
expended  in  firing  the  I^1  mission*. 

Parameters  (9)  through  (13)  are  all  computed  by 
subroutine  SERVIC,  which  will  be  discussed  in  Section  III. 6. 

(14)  FSTSHT(I) :  The  time  between  the  call  for 
fire  and  the  time  the  first  round  was  fired  for  a  mission 
which  required  no  queue  time.  FSTSHT(I)  *  TRAVEL(I)  +  SETUP(I). 


(15)  FSTSHO(I) :  Same  as  FSTSHT(I)  for  a  mission 
which  required  queue  time.  FSTSHO(I)  *  QTIMB(x)  ♦  TRAVEL f I) 

♦  SETUP (I). 

(16)  JOINQ(I) :  The  time  the  Ith  mission  joined 
the  queue.  JOINQ(I)  s  0  if  the  I**1  mission  required  no  queue 
time . 

(17)  LEAVEQ(I):  Time  the  Ith  mission  left  the 
queue. LEAVEQ(I)  *  0  if  the  Ith  mission  required  no  queue  time. 

(18)  QTIME(I):  Amount  of  time  the  Ith  mission 
was  required  to  wait  for  assignment  of  a  ship  to  fire  the 
mission  (queue  time).  QTIME(I)  »  LEAVEQ(I)  -  JOINQ(I). 

(19)  MSNCMP(I) :  The  time  at  which  the  Ith 
mission  was  completed.  MSNCMP(I)  *  MSNSKD(I)  +  QTIME(I) 

*  SERVC(I) . 

D.  SHIP  ASSIGNMENT  DOCTRINE 

1.  Types  of  Ship  Assignment  Doctrines 

There  are  two  types  of  ship  assignment  doctrine 
available  to  the  user.  These  refer  to  the  decision  process 
which  the  STDF  commander  uses  to  determine  which  types  of 
ships  from  among  those  "present  in  the  problem"  are  appro¬ 
priate  to  fire  a  particular  mission.  The  type  of  ship  as¬ 
signment  doctrine  is  specified  by  the  input  variable  ASNDOC , 
which  has  values  defined  as  follows: 

1 :  "Any  ship  with  guns  of  adequate  range  may  fire. " 
The  target  range  is  the  only  consideration. in  determining 
which  types  of  ships  are  appropriate  to  assign  to  fire  the 
mission. 
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2:  "Best  effect  on  target  ♦  "  Both  target  range 
and  target  type  are  considered  in  determining  which  types 
of  ships  are  appropriate  to  assign  to  fire  the  mission. 

When  ASNDOC  *  2  is  used,  the  following  rules  govern 
the  determination  of  appropriate  ships: 

^  1.  Heavy  material  targets  (type  3)  may  only  be 

fired  at  by  8"/55  or  6"/47  ships,  unless  there  are  no  8"/55 
or  6"/47  ships  present  in  the  problem  or  there  are  no  8"/5S 
ships  present  and  the  target  is  out  of  range  of  a  6"/47  gun. 
The  reason  for  this  rule  is  that  8"/55  and  6"/47  shells  have 
a  considerably  higher  explosive  power  than  5"/54  or  S"/38 
shells. 

2.  Personnel  targets  (types  1  or  2)  may  only  be 
fired  at  by  5"/54  or  S'738  ships  and  a  5"/54  ship  is  prefer- 
rable  due  to  its  rapid  rate  of  fire,  unless  there  are  no 
5"/54  or  5*738  ships  present,  or  the  target  is  out  of  range 
of  a  5"/54  gun,  or  there  are  no  5"/54  ships  present  and  the 
target  is  out  of  range  of  a  5"/38  gun. 

2.  Mission  Types 

Regardless  of  which  assignment  doctrine  is  used, 
each  mission  is  assigned  a  mission  type  number  (MISTYP(I)). 
This  number  indicates  what  types  of  ships,  from  among  those 
present  in  the  problem,  are  appropriate  (eligible)  to  fire 
the  mission.  There  are  8  possible  values  of  the  variable 
MISTYP(I) ,  as  follows: 


41 


MTCTVD  f"M 

1 

2 

3 

4 

5 

6 

7 

8 


Twno  Cli  4 «b  Cl  4  #»4Ul  a  Cl 


8‘7SS 

8'*/S5,  6*'/47 
5*'/54 

5’’/54,  5**/38 
Any 

8’75S,  6H/47,  5"/54 
8"/55 ,  5M/54 
None 


3.  Maximum  Gun  Ranges 

The  maximum  range  of  each  type  gun  in  the  problem 
is  specified  by  the  user  as  the  input  variables  MAX855, 
MAX647,  MAXS54,  and  MAX538.  The  maximum  range  of  the  longest 
type  of  gun  present  is  also  specified  by  the  input  variable 
MAXNGF.  Due  to  the  construction  of  the  model  there  are  two 
requirements  on  these  quantities: 

1.  MAX855  <  MAX554  <  MAX647  <  MAX538,  unless  no 
ships  of  a  particular  type  are  present  (i.e.,  if  NR554  *  0, 
then  MAX554  =  0) . 

2.  MAXNGF  <  MAX  {NORTH (NR) > 

NR 

In  view  of  the  first  requirement,  the  only  type  of 
missions  which  will  occur  if  ASNDOC  *  1  are  1,5,6,  and  7. 

Any  of  the  types  1  through  7  may  occur  if  ASNDOC  *  2.  A 
type  8  mission  will  never  occur  if  the  second  requirement 
is  satisfied.  If  a  type  8  mission  does  occur,  the  message 
**A  TYPE  8  MISSION  HAS  OCCURRED  IN  AREA  will  be  printed 

out  and  the  program  will  terminate. 
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4.  Warning  Concerning  Use  of  Ship  Assignment  Doctrine  2 
Assignment  doctrine  2  has  one  definite  shortcoming. 
It  only  considers  which  types  of  ships  are  present  in  the 
problem,  not  which  types  are  available  to  fire  at  the  time 
the  mission  is  called  for.  Thus,  a  situation  could  occur 
in  which  the  only  types  of  ships  appropriate  to  firr.  a 
mission  are  off-station  for  rearming  or  refueling  at  the 
time  the  mission  is  call'ed  for.  Then  the  mission  would  be 
required  to  wait  until  one  of  these  ships  returns  before  it 
is  fired  although  other  ships  with  guns  of  adequate  range 
might  be  available  probably,  resulting  in  an  unacceptably 
large  amount  of  queue  time. 

E.  SHIP  AVAILABILITY 

After  a  mission  has  been  assigned  a  mission  type  in  ac¬ 
cordance  with  the  rules  of  the  specified  assignment  doctrine, 
subroutine  AVAIL  is  called  once  for  each  type  of  ship  which 
is  eligible  to  fire  the  mi ssion  and  present  in  the  problem. 
The  arguments  of  AVAIL  are  the  number  of  ships  present  of 
the  type  being  considered,  and  the  type  of  ship  being  con¬ 
sidered. 

1.  Ship  Status 

The  status  of  each  ship  in  the  problem  is  maintained 
by  the  variable  STATUS (S,T),  where  S  is  the  ship  number  and 
T  is  the  ship  type.  There  are  four  possible  values  of 
STATUS (S,T),  as  follows: 


1:  Ship  available  for  assignment. 

2:  Ship  presently  assigned  to  a  mission. 


3:  Ship  presently  refuel'  or  in  a  "refuel  after 
mission"  status,  or  engaged  in  a  combined  refueling/rearming. 

4:  Ship  presently  rearming. 

2 .  Vector  of  Available  Ships 

The  series  cf  calls  on  AVAIL  for  each  mission  re¬ 
sults  in  a  vector  (/AIL)  of  ships  of  types  appropriate  for 
the  mission  for  which  STATUS (S,T)  *  1.  The  values  in  this 
vector  are  two-digit  numbers  with  the  ship  number  as  the 
first  digit  and  ship  type  as  the  second  digit.  For  example, 
VAIL(l)  =  23  means  that  the  first  ship  checked  which  was 
available  was  5"/54  ship  #2.  S/R  AVAIL  also  keeps  a  tally 
of  the  number  of  ships  which  are  available  to  fire  the  mis¬ 
sion.  This  tally  is  the  value  of  the  variable  KOUNT. 

3.  Assignment  of  Closest  Ship 

•-’2r  the  vector  of  available  ships  is  obtained  by 
call'  ■'*  AVAIL,  control  is  returned  to  the  main  program.  If 
KOUNT  =  0,  the  mission  joins  the  queue.  If  KOUNT  >  0,  then 
the  distance  of  each  ship  in  the  vector  VAIL  from  the  firing 
position  is  determined  in  the  following  manner: 

Each  available  sh:p  is  considered  to  be  at  some 
random  point  in  the  waiting  pattern.  For  each  ship,  its 
distance  from  the  firing  position  is  calculated  as  a  uni¬ 
formly  distributed  random  variable  between  the  m.  nimum  pos¬ 
sible  distance  (MINDIS)  and  the  maximum  possible  distance 
(M/XDIS) .  The  method  of  calculation  of  MINDIS  and  MAXDIS 
is  indicate1  in  Figure  2. 

The  closest  ship  is  then  assigned  to  fire  the  mis¬ 
sion.  Thi:-  is  done  by  the  following  calculations: 
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MAXDIS  -  VfWTDIS+NSWAIT) 2  +  [MAX (WTNE - EWTAR ( I) , EWTAR C I ) -  WTNW) ] 
Figure  2.  Method  of  Calculation  of  MINDIS  and  MAXDIS. 


TYPSHP(I)  *  MOD  (VAIL(JD) ,  10) 

HULLNR(I)  =  VAIL (JD) /IQ 

*  1 

where  JD  is  the  position  in  the  vector  VAIL  of  the  closest 
ship. 

F.  THE  QUEUE 

The  queue  in  the  model  consists  of  those  missions  for 
which  the  time  of  call  of  fire  is  less  than  the  current  prob¬ 
lem  time  (MSNSKD(I)  <  CLOCK)  but  to  which  a  ship  has  not  yet 
been  assigned.  Thus,  missions  in  the  queue  are  referred  to 
as  ’'awaiting  assignment”. 

1 .  Joining  the  Queue 

A  mission  joins  the  queue  when  no  ship  appropriate 
for  the  mission  type  is  available  at  the  time  the  mission  is 
called  for  (KOUNT  =  0) .  Mission  numbers  of  missions  awaiting 
assignment  are  maintained  in  the  vector  Q.  The  current  num¬ 
ber  of  missions  awaiting  assignment  is  maintained  by  the  vari 
able  INQ.  In  view  of  the  assumption  stated  in  paragraph  II. 
D.S,  the  maximum  possible  number  of  missions  in  the  queue  is 
equal  to  the  number  of  patrol  subareas. 

2 .  Leaving  the  Queue 

A  mission  leaves  the  queue  when  an  appropriate  type 
ship  becomes  available  by  having  completed  another  mission  or 
having  completed  refueling  and/or  rearming,  and  there  are  no 
other  missions  ahead  of  it  in  the  queue  which  are  waiting 
for  the  same  type  of  ship.  Each  time  a  ship  becomes  avail¬ 
able  and  INQ  >  0,  each  mission  in  the  queue,  beginning  with 
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Q(l)  is  checked  to  see  if  the  available  ship  is  appropriate 
to  that  mission.  As  soon  as  one  is  located,  that  mission  is 
removed  from  the  queue  and  all  other  missions  following  it 
are  moved  up  one  number  in  the  queue. 

3.  Queue  Discipline 

There  are  two  types  of  queue  discipline  available 
to  the  user.  The  type  desired  is  specified  by  the  input 
variable  QDISC,  which  may  have  the  following  values: 

1:  First-come  first  served. 

2:  Personnel  missions  have  priority. 

The  use  of  QDISC  =  2  places  type  1  and  2  missions  ahead  of 
type  3  and  4  missions  in  the  queue.  There  is  no  priority 
relationship  between  type  1  and  2  missions.  Nc  e  that 
QDISC  »  2  should  not  be  used  when  ASNDOC  =  2,  or  vice-versa, 
since  ASNDOC  =  2  only  allows  5"/38  ships  to  fire  personnel 
missions.  There  is  no  provision  in  the  program  for  missions 
which  are  already  in  progress  (i.e.,  a  ship  has  been  assigned) 
to  be  preempted  by  any  other  mission.  Thus,  once  a  ship  is 
assigned  to  fire  a  mission,  it  is  not  available  for  any  other 
mission  until  that  mission  is  completed. 

G.  SERVICE  TIME 

Each  time  a  ship  is  assigned  to  fire  a  mission,  subrou¬ 
tine  SERVIC  is  called  to  calculate  the  service  time  for  that 
mission.  The  arguments  of  SERVIC  are  type  of  ship  assigned, 
target  type,  target  range,  mission  number,  and  distance  of 
assigned  ship  from  the  firing  position.  S/R  SERVIC  calcu¬ 
lates  the  total  service  time  of  the  mission  and  amount  of 
ammunition  expended  during  the  mission. 
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1.  Total  Service  Time  (SERVCCD) 

There  are  three  components  of  total  service  time: 
a.  Travel  Time  (TRAVEL (I)) 

Travel  time  is  the  time  required  for  the  ship 
to  reach  the  firing  position  after  it  has  been  assigned.  It 
is  calculated  in  one  of  three  ways. 

(1)  Ship  Coming  from  the  Waiting  ^lttern.  If 
the  assigned  ship  was  assigned  by  virtue  of  having  been  the 
closest  available  ship  of  appropriate  type  at  the  time  the 
mission  was  called  for,  then  its  travel  distance  is  the  dis¬ 
tance  calculated  by  the  method  described  in  Section  III.E.3. 

(2)  Ship  Having  Just  Completed  Another  Mission. 
In  this  case  the  travel  distance  is  the  distance  along  the 
gunline  between  the  two  firing  positions.  This  distance  is 

RUN  *  IABS  (EWTAR(E) -EWTAR(P) ) 

where  E  is  the  number  of  the  mission  just  completed  and  P  is 
the  number  of  the  new  mission  to  which  the  ship  has  been  as¬ 
signed. 

(3)  Ship  Having  Just  Completed  Refueling  and/or 
Rearming.  A  ship  is  considered  to  have  completed  refueling 
and/or  rearming  when  it  returns  to ^a.  position  midway  between 
the  northern  and  southern  edges  of  the  waiting  pattern  and 
due  south  of  the  firing  position.  Then  travel  distance  is 

RUN  =  WTDIS  +  NSWAIT/2 . 

Regardless  of  which  of  the  above  methods  is  used 
to  calculate  travel  distance,  S/R  SERVIC  calculates  travel 
time  as 
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TRAVEL (I)  =  DIS/ (SPEED  *  33.33) 


where  DIS  is  the  appropriate  travel  distance,  SPEED  is  an 
input  variable  which  specifies  the  speed  (in  knots)  at  which 
ships  proceed  to  the  firing  position,  and  33.33  is  a  factor 
which  converts  knots  (nautical  miles/hour)  to  yards/minute . 

b.  Setup  Time  (SETUP (I)) 

This  is  the  time  required  for  a  ship  to  fix  its 
navigational  position  and  ready  its  guns  after  it  has  reached 
the  firing  position.  S/R  SERVIC  calculates  setup  time  as  a 
uniformly  distributed  random  variable  between  the  minimum  ex¬ 
pected  setup  time  (MINSET)  and  maximum  expected  setup  time 
(MAXSET) ,  where  MINSET  and  MAXSET  are  input  variables  spec¬ 
ified  by  the  user. 

c.  Firing  Time  (FIRE(I)) 

Actual  firing  time  is  the  time  interval  between 
the  time  the  first  round  is  fired  and  the  time  the  last  round 
is  fired.  For  each  mission  this  is  considered  to  be  an  ex¬ 
ponentially  distributed  random  variable  with  a  mean  which  is 
a  function  of  the  type  of  ship  firing,  the  target  type,  and 
the  target  range. 

(1)  Mean  Firing  Time  Due  to  Target  Range.  The 
portion  of  the  mean  firing  time  which  is  due  to  the  target 
range  is  considered  to  be  a  linear  function  of  the  target 
range.  This  portion  is  calculated' as  RANGE(I)/DENOM(J,K) , 
where  J  is  the  type  ship  and  K  is  the  target  type.  The 
values  of  DENOM(J,K)  are  specified  by  the  user  for  each  pos¬ 
sible  combination  of  ship  type  vs.  target  type.  Thus,  DENOM 
is  a  4  x  4  matrix. 
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(2)  Mean  Firing  Time  Due  to  Target  Type.  This 
quantity  is  specified  by  the  user  as  FIRTIM(J,K)  where  J  is 
the  ship  type  and  K  is  the  target  type.  Thus,  FIRTIM  is 
also  a  4  x  4  matrix. 

The  mean  of  the  exponentially  distributed  firing 
time  for  each  mission  is  calculated  as 

MU  =  RANGE ( I ) / DENOM ( J , K)  +  FIRTIM(J,K). 

Therefore,  the  firing  time  for  each  mission  is  calculated  as 

FIRE(I)  =  -MU*ALOG(XRN) . 

After  travel,  setup,  and  firing  time  have  been 
calculated,  total  service  time  is  calculated  as 

SERVC(I)  =  TRAVEL ( I )  +  SETUP(I)  +  FIRE(I). 

2.  Ammunition  Expenditure  (AMMEXP(I)) 

S/R  SERVIC  also  calculates  the  amount  of  ammunition 
expended  during  each  mission.  The  user  specifies  the  minimum 
(MINAMM(J,K))  and  maximum  (MAXAMM ( J , K) )  expected  rounds  of 
ammunition  expended  per  minute  by  ship  type  J  firing  against 

target  type  K.  These  values  should  take  into  account  the 

'  y 

average  rounds  per  minute  over  an  entire  mission,  i.e.,  the 
slow  rate  of  fire  at  the  beginning  of  the  mission  and  the 
rapid  rate  of  fire  after  the  fire  has  been  spotted  onto  the 
target  should  be  averaged.  Therefore,  MAXAMM (J,K)  should 
normally  not  be  the  maximum  attainable  rate  of  fire  by  a 
type  J  ship. 

The  number  of  rounds  of  ammunition  expended  per  minute 
for  each  mission  is  calculated  as  a  uniformly  distributed 
random  variable  between  MINAMM(J,K)  and  MAXAMM(J,K).  Therefore 
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the  number  of  rounds  of  ammunition  expended  during  the  I*h 
mission  is  calculated  as 

AMMEXP(I)  =  (MINAMM(J,K)+XRN*(MAXAMM(J,K)-MINAMM(J,K)))*FIRE(I) 
where  J  is  the  type  of  ship  and  K  is  the  type  of  target. 

H.  REFUELING  AND  REARMING 

The  user  has  two  options  available  for  taking  refueling 
and  rearming  of  ships  into  account  in  the  problem.  The  op¬ 
tion  desired  is  specified  by  the  input  variable  REPPOL  (Re¬ 
plenishment  Policy),  which  has  values  defined  as  follows: 

1:  Ships  do  not  refuel  or  rearm. 

2:  Ships  refuel  according  to  a  predetermined  schedule 
and  rearm  when  their  on-board  ammunition  level  reaches  a 
specified  minimum  level. 

The  case  in  which  ships  do  not  rearm  't>r  refuel  would 
arise  when  enough  ships  are  available  to  replace  ships  which 
are  required  to  depart  for  refueling  and/or  rearming,  thus 
the  departure  of  the  affected  ship  would  not  effect  the 
queueing  results. 

If  REPPOL  =  2,  refueling  and  rearming  is  handled  in  the 
following  ways  in  the  program: 

1.  Refueling  * 

a.  Methods  of  Establishing  a  Refueling  Schedule 
The  user  establishes  a  refueling  schedule  for 
all  ships  for  the  entire  duration  of  the  problem.  There  are 
two  options  available  for  establishing  this  schedule.  The 
desired  option  is  specified  by  the  input  variable  REFSKD, 
which  has  values  defined  as  follows: 
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1:  The  user  specifies  the  exact  time  of  each 
scViduied  refueling  for  each  ship.  For  each  type  of  ship,  a 
data  card  is  prepared  which  contains  the  floating  point  time 
(in  days  from  zero  problem  time  to  2  decimal  places)  that 
each  ship  of  that  type  is  to  depart  station  for  its  first 
scheduled  refueling.  Similar  cards  are  prepared  for  the 
2nd,  3rd,  and  up  to  the  9th  scheduled  refuelings  for  ships 
of  each  type.  These  refueling  times  are  the  values  of  the 
variable  XFULSKfl ,J,X) ,  where  I  is  the  ship  number,  J  is 
the  ship  type,  and  K  is  the  number  of  the  refueling.  The 
total  number  of  rounds  of  refueling  is  specified  by  the  in¬ 
put  variable  FINFUL  (must  be  <9) .  Thus,  there  must  be  4  x 
FINFUL  cards  in  the  refueling  schedule  data  deck.  The  pro - 
gram  converts  these  scheduled  refueling  times  to  minutes 
from  zero  problem  time  and  places  them  on  the  Event  Chain  in 
the  proper  chronological  order. 

2:  The  refueling  schedule  scheme  discussed 
above  is  replaced  by  a  simpler  (although  less  flexible) 
scheme.  The  entire  schedule  is  specified  on  one  data  card. 
The  first  five  quantities  on  the  card  are  the  floating  point 
times  (in  days  from  zero  problem  time  to  2  decimal  places) 
that  the  first  cruiser  is  to  depart  for  its  1st  thru  Sth 
scheduled  refueling.  These  are  specified  by  the  input  vari¬ 
ables  XCRFST(I) ,  1=1  thru  5.  The  next  5  quantities  are  the 
times  the  first  destroyer  is  to  depart  for  its  1st  through 
5th  scheduled  refuelings.  These  are  specified  by  the  input 
variables  XDDFST(I),  1=1  thru  5.  The  last  two  numbers  are 
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the  time  interval  in  minutes  at  which  succeeding  cruisers 
and  destroyers  follow  the  first  for  refueling.  These  are 
specified  by  the  input  variables  CRINT  and  DDINT.  The  pro¬ 
gram  then  calculates  a  refueling  schedule  for  the  entire 
duration  of  the  problem  using  these  parameters.  There  are 
several  rules  which  must  be  followed  in  the  preparation  of 
the  refueling  schedule  data  cards. 

a.  For  ships  which  are  not  present  in  the 
problem,  the  appropriate  columns  on  the  data  cards  should  be 
left  blank. 

b.  If  REFSKD  *  1,  for  ships  which  are  present 
in  the  problem  but  are  not  scheduled  to  refuel  FINFUL  times, 
a  large  number  (99.00  for  example)  must  be  inserted  in  the 
appropriate  columns  on  the  data  cards. 

c.  If  REFSKD  =  2  and  less  than  five  rounds  of 
refueling  are  to  be  used,  a  large  number  (99.00)  must  be  in¬ 
serted  in  the  appropriate  columns  on  the  data  card. 

d.  If  REFSKD  =  1  and  no  ships  of  a  particular 
type  are  present  in  the  problem,  a  blank  card  must  be  in¬ 
serted  in  the  data  deck  in  the  appropriate  place  for  each 
round  of  refuelings. 

Examples 

REFSKD  =  1 :  The  expected  duration  of  the  problem  is  10.0 
days.  NR855  =  2,  NR647  =  0,  NR554  =  1,  NR538  =2.  It  is 
desired  that  the  cruisers  refuel  on  the  sixth  day  and  the 
destroyers  refuel  on  the  fifth  and  ninth  days.  Then  typical 
values  of  XFULSK(I,J,K)  would  be 
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First  Card 

Fifth  Card 

XFULSK (1,1,1) 

=  5.50 

XFULSK(1 ,1 ,2) 

=  99.00 

XFULSK(2 ,1,1) 

=  5.75 

XFULSK(2 ,1 ,2) 

=  99.00 

Third  Card 

Seventh  Card 

XFULSK(1,3,1) 

=  4.25 

XFULSK (1,3, 2) 

=  8.25 

Fourth  Card 

Eighth  Card 

XFULSK(1 ,4,1) 

=  4.50 

XFULSK(1 ,4 ,2) 

=  8.50 

XFULSK (2, 4,1) 

=  4.75 

XFULSK (2, 4, 2) 

*  8.75 

FINFUL  =  2(2  rounds  of  refueling) .  Therefore  there  would 
be  8  data  cards,  of  which  the  2nd  and  6th  data  cards  would 
be  blank  since  NR647  =  0. 

REFSKD  =  2 :  For  the  same  situation  described  above,  the  fol- 
lowing  values  on  one  data  card  would  produce  the  same  schedule 
as  above: 

XCRFST (1)  =5.50  XDDFST(l)  =4.25 

XCRFST (2)  thru  XCRFST(5)  =  99.00  XDDFST(2)  =  8.25 

CRINT  =360  *  XDDFST(3)  thru 

XDDFST(5)  =  99.00 

DDINT  =  360 

b.  Latest  Fueling  Time 

The  program  calculates  the  time  of  extraction 

of  each  patrol  from  its  subarea  as  INSEIT(NR)  plus  a  uni- 

*  ' 

formly  distributed  random  variable  between  XSHORT(NR)  and 
XLONG(NR) .  The  latest  of  these  extraction  times  is  stored 
as  the  variable  MAXTRC.  In  computing  the  fueling  schedule 
according  to  user  input  parameters,  no  refuelings  are  placed 
on  the  Event  Chain  which  have  a  scheduled  time  later  than 
MAXREF  =  MAXTRC  -  NOMO,  where  NOMO  is  an  input  variable 
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specified  by  the  user.  MAXREF  is  defined  as  the  latest  time 
that  ships  will  depart  for  refueling.  For  example,  if  it 
is  desired  that  no  ship  depart  for  refueling  la*er  than  6 
hours  before  the  latest  anticipated  extraction  time,  then 
NOMO  *  360. 

2.  Rearming 

The  user  specifies  several  input  variables  which 
are  used  to  determine  when  each  ship  requires  rearming.  These 
input  variables  are  defined  as  follows: 

INAMM(S,T):  The  initial  ammunition  level  (rounds 
on  board)  of  ship  number  S  of  type  T. 

MINAMM(J,K)  and  MAXAMM(J,K) :  Defined  in  Section 

III. 6.2. 

MINLEV (T) :  The  minimum  ammunition  level  to  which  a 
ship  of  type  T  is  allowed  to  fall  before  it  is  required  to 
rearm. 

Each  time  a  mission  is  fired,  ammunition  expenditure 
for  that  mission  is  calculated  by  S/R  SERVIC  as  discussed  in 
Section  III.G.2.  This  amount  is  subtracted  from  the  current 
ammunition  level  of  the  ship  vhich  fired  the  mission,  as  fol¬ 
lows  : 

AMMLEV(S,T)  *  AMMLEV (S ,T)  -  AMMEXP(I). 

If  AMMLEV(S,T)  <  MINLEV (T) ,  the  ship  is  sent  to  rearm. 

3.  Rearming  and  Refueling  Completions 

a.  Refueling 

When  an  event  for  which  EVTYP(K)  =  3  (scheduled 
refueling)  occurs,  the  status  of  the  scheduled  ship  is  checked. 
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If  the  ship  is  presently  firing  a  mission  (STATUS (S,T)  *  2) 
it  is  placed  in  a  “refuel  after  mission*'  status  (STATUS(S,T) 

13  3) .  If  STATUS(S,T)  =  1,  the  ship  is  sent  to  refuel  im¬ 
mediately.  In  either  case,  the  time  of  completion  of  fueling 
is  calculated  as  follows: 

The  user  specifies  the  minimum  and  maximum  ex¬ 
pected  lengths  of  time  that  it  should  take  for  a  cruiser  and 
a  destroyer  to  refuel.  These  are  the  input  variables  MINFLC, 
MAXFLC,  MINFLD,  and  MAXFLD.  These  times  should  take  into 
account  travel  time  to  and  from  the  rendezvous  point  as  well 
as  time  alongside  the  oiler.  For  example,  if  it  is  planned 
that  ships  will  rendezvous  with  the  oiler  at  a  point  20 
miles  from  the  center  of  the  waiting  pattern  and  will  proceed 
to  the  rendezvous  point  at  a  speed  of  20  knots,  and  that  a 
destroyer  should  require  30  to  60  minutes  alongside  the 
oiler,  than  MINFLD  should  be  60+30+60  *  150,  and  MAXFLD 
should  be  60+60+60  =  180. 

The  time  of  completion  of  refueling  (COMFUL)  is 
calculated  as  CLOCK  time  (the  time  the  ship  departed  for  re¬ 
fueling,  either  after  a  mission  is  completed  or  at  the  sched¬ 
uled  refueling  time)  plus  a  uniformly  distributed  random 
variable  between  MINFLC  and  MAXFLC  for  cruisers,  or  MINFLD 
and  MAXFLD  for  destroyers.  COMFUL  is  then  placed  on  the 
Event  Chain. 

b.  Rearming 

Completion  of  rearming  (COMARM)  is  handled  in 
the  same  way  as  COMFUL  except  that  the  parameters  of  the 
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uniform  distribution  of  rearming  time  are  MINARC,  MAXARC, 

MI  'ARD,  and  MAXARD,  specified  by  the  user  with  the  same 
considerations  as  the  refueling  time  parameters.  When 
COMARM  is  calculated,  the  rearming  ship's  ammunition  level 
(AMMLEV(S,T))  is  set  equal  to  its  initial  ammunition  level 
(INAMM(S,T) ) . 

c.  Combined  Re fueling/ Re arming 

When  a  ship  which  is  scheduled  to  refuel  is 
off-station  for  rearming  at  the  time  the  scheduled  refueling 
occurs,  then  a  refueling  time  calculated  in  the  manner  de¬ 
scribed  above  is  added  to  the  previously  calculated  COMARM. 
From  this  time,  the  input  variable  COMBO  is  subtracted. 

COMBO  is  defined  as  the  amount  of  time  off-station  saved  by 
conducting  both  refueling  and  rearming  in  the  same  evolu¬ 
tion.  Thus,  it  should  be  approximately  equal  to  the  time 
for  a  one-way  trip  from  the  ammunition  ship  to  the  waiting 
pattern  plus  the  time  for  a  one  way  trip  from  the  waiting 
pattern  to  the  oiler,  minus  the  time  required  to  travel  from 
the  ammunition  ship  to  the  oiler.  COMBO  must  not  be  greater 
than  MINFLD  or  MINFLC . 

The  previously  calculated  COMARM  is  removed 
from  the  Event  Chain  and  the  new  time  of  completion  of  com¬ 
bined  refueling/rearming  (COMREP)  is  placed  on  the  Event 
Chain. 

I .  OUTPUT  DATA 

All  output  of  the  program  is  well-formatted  to  facili¬ 
tate  reading  and  interpretation.  Each  type  of  output  is  dis¬ 
cussed  below. 
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1 .  Printout  of  Input  Data 

A  printout  of  the  user- specified  values  of  most  of 
the  input  variables  is  produced  at  the  beginning  of  the  out¬ 
put  data.  Appendix  C  contains  an  example  of  this  printout. 

2 .  Mission  History 

The  Mission  History  is  a  chronological  list  of  the 
12  most  significant  parameters  assigned  to  each  mission  dur¬ 
ing  a  particular  iteration  of  the  problem.  There  are  three 
options  available  to  the  user  concerning  the  Mission  History 
printout.  The  desired  option  is  specified  by  the  input 
variable  MISHIS,  which  has  values  as  follows: 

1:  No  Mission  History  printout  desired. 

2:  Printout  Mission  History  for  the  first  iteration 

only. 

3:  Printout  Mission  History  for  each  iteration. 

The  time  of  call  for  fire  and  completion  of  each  mission  are 
converted  to  date-time  groups  in  the  Mission  History  printout. 
The  first  1  or  2  digits  are  the  problem  day  and  the  last  four 
digits  are  the  24 -hour  clock  time. 

Appendix  D  contains  an  example  of  the  Mission  His¬ 
tory  printout. 

3.  Statistical  Results  for  Each  Iteration 

For  each  iteration,  a  printout  of  the  pertinent 
queueing,  ship  utilization,  refueling,  rearming,  and  ammu¬ 
nition  expenditure  statistics  is  available.  These  are  the 
results  for  that  iteration  only,  they  are  not  cumulative 
over  all  iterations  up  to  that  point.  The  user  specifies 
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whether  printout  of  iteration  statistical  results  is  desired 
by  the  input  variable  ITSTAT,  which  may  have  the  following 
values : 

1  •  Drintout  of  iteration  statistical  results  is 

desired. 

2  **  Printout  of  statistical  results  for  each  itera¬ 
tion  is  not  desired. 

Appendix  E  contains  an  example  of  the  printout  of 
iteration  statistical  result;. 

4.  Overall  Statistical  Results 

After  all  NRIT  iterations  of  the  problem  have  been 
completed,  a  printout  of  the  pertinent  overall  statistical 
results  is  produced.  This  printout  includes  the  mean  and 
the  standard  deviation  of  the  mean  for  each  quantity  except 
for  refueling,  rearming,  and  ammunition  expenditure  data, 
for  which  the  mean  only  is  printed  out. 

In  computing  the  mean  and  the  standard  deviation 
of  the  mean  of  the  quantities  in  the  overall  statistical  re¬ 
sults,  the  following  formulas  are  used: 


1  n 

wheie  is  the  result  of  the  it}l  iteration  and  n  =  NRIT. 

Using  the  means  and  standard  deviation  from  the 
overall  statistical  results  for  NRIT  iterations,  confidence 
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intervals  tor  the  actual  value  of  each  output  quantity  can 
be  calculated  using  the  following  formula: 


0x 


t(n-l) 

(l-cO/2 


<  M  <  X  +  0= 


t(n-l) 

tCi-a)/2 


where  x  is  the  observed  mean,  a=  is  the  observed  standard  de- 
viation  of  the  mean,  and  *s  t^ie  value  from  a  table 

values  of  the  Student 's-t  distribution  with  n-1  degrees 
of  freedom  at  the  (l-ct)/2  level  of  significance .  a  is  de¬ 
fined  as  the  probability  that  the  actual  value  of  the  quan¬ 
tity  lies  outside  of  a  confidence  interval  calculated  by 
the  above  f  rmula.  Therefore  (1-a)  is  the  degree  of  con¬ 
fidence  with  which  statements  concerning  the  actual  value 
of  the  quantity  can  be  made. 

Example 

U  =  AVERAGE  WAITING  TIME  PER  MISSION 

10  iterations  of  a  problem  yield  x  =  20.0  and  a=  =  1.00.  The 
user  desires  to  be  able  to  predict  with  95$  confidence  the 
interval  in  which  the  actual  value  of  AVERAGE  WAITING  TIME 
PER  MISSION  lies  based  on  these  results  of  10  iterations. 

A  table  of  values  of  the  Student’ s-t  distribution  shows 
=  2.26.  Therefore,  it  can  be  stated  with  95$  con¬ 
fidence  that  1 


20.0  -(1.00) (2.26)  <  p  <  20.0  +(1.00) (2,36) 


or 


17.74  <  AVERAGE  WAITING  TIME  PER  MISSION  <  22.26. 
Appendix  F  contains  an  example  of  the  printout  of 
overall  statistical  results. 
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IMPIfT  DATA  ~  niTi  ncr*/ 

—  “ '  “  —  -  — **  •*  “■*  ■*  *  x-v  t/uvn, 


This  section  lists  all  input  variables  in  the  data  deck 
by  card  and  column  number(s) .  For  each  variable  the  follow¬ 
ing  information  is  listed:  variable  name,  subsection  of 
section  III  in  which  the  variable  is  described,  the  FORTRAN 
variable  type,  the  maximum  possible  value  which  can  be  as¬ 
signed,  and  the  units  of  the  indicated  quantity.  All  quanti¬ 
ties  must  be  right-adjusted  within  their  appropriate  columns. 


1 ‘  |i*st  Card:  Miscellaneous  Problem  Definition  Vari- 


Cols 

Variable 

Re ference 

Type 

Maximum 

Units 

1-2 

NHAREA 

B.l 

12 

50 

Subareas 

3 

NR855 

A.  S 

11 

9 

Ships 

4 

NR647 

A.  5 

11 

9 

Ships 

5 

NR554 

A.  5 

11 

9 

Ships 

6 

NR538 

A. 5 

11 

9 

Ships 

7-11 

MAX8S5 

D.3 

IS 

99999 

Yords 

12-16 

MAX647 

D.3 

15 

99999 

Yards 

17-21 

MAXS54 

D.  3 

IE 

99999 

Yards 

22-26 

MAXS38 

D.3 

15 

99999 

Yards 

27-31 

MAXNGF 

D.3 

15 

99999 

Yards 

32-38 

WTNW 

B.3 

17 

9999999 

Yards 

39-45 

WTNE 

B.3 

17 

9999999 

Yards 

46-51 

NS WAIT 

B.3 

16 

999999 

Yards 

52-57 

WTD  IS 

B.3 

16 

999999 

Yards 

58-59 

SPEED 

G.  1 .  a 

12 

99 

Knots 

60-61 

MINSET 

G.l.b 

12 

99 

Minutes 

62-63 

MAXSET 

G.l.b 

12 

99 

Minutes 

64-67 

HHOUR 

C.2.a. (1) 

14 

1439 

Minutes 
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2. 

=■*■=»  = 

vvWHa  =_-a. 

rd’  Program 

Control 

In4i?a*or 

Cpli 

Variable 

Re  feT/jncp 

Type 

Xtn  v  4  mt*M 

*  **“*>  *  fllWHi 

Hn  U. 

WH  IV? 

1 

ASNDOC 

D.l 

11 

2 

2 

QDISC 

F.  3 

11 

2 

— 

3 

REPPOL 

H 

11 

2 

— 

4 

REFSKD 

H.l 

11 

2 

— 

5 

FINFUL 

H.l 

11 

9 

Rounds  of 
Refueling 

6 

MI SHIS 

1.2 

11 

3 

— 

7 

ITSTAT 

1.3 

11 

2 

8-10 

NRIT 

A. 7 

13 

999 

Iterations 

3*  Next  NRAREA  Cards;  Subarea  Information  Cards 


One  card  is  prepared  for  each  patrol  subarea.  Cards 


are  placed  in  the  data  deck  in 

is  the  subarea  number. 

order  of 

subarea 

numbers .  NR 

Cols 

Variable 

Reference 

im. 

Maximum 

Units 

1-7 

NORTH (NR) 

B.2 

17 

9999999 

Yards 

8-14 

SOUTH (NR) 

B.2 

17 

9999999 

Yards 

15-21 

EAST (NR) 

B.l 

17 

9999999 

Yards 

22-28 

WEST (NR) 

B.l 

17 

9999999 

Yards 

29-33 

XCFF(NR) 

C.  2 .  a 

F5.1 

999.9 

Missions 

34-38 

XSHORT(NR) 

C.2.a 

F5.1 

999.9 

Days 

39-43 

XLONG(NR) 

C.2 .  a 

F5.1 

999.9 

Days 

44-49 

INSERT (NR) 

C.2 .a. (1) 

16 

999999 

Minutes 

50-54 

XPROB1  (NR) 

C .  2 .  b .  (  5  ) 

F5.2 

1.00 

Probability 

55-59 

XPR0B2 (NR) 

C.2 .b. (5) 

F5.2 

1.00 

Probability 

60-64 

XPR0B3(NR) 

C.2.b. (5) 

F5.2 

1.00 

Probability 

65-69 

XPR0B4 (NR) 

C.2.b.  (5) 

F5.2 

1.00 

Probability 
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First 

Card:  Mean  Firing  Time 

T\l*  A  ^  A 
DVC  UU 

l  a*j^v  w 

Cols 

Variable 

Re  ference 

Type 

Maximum 

Units 

1-2 

FIRTIM(1 ,1) 

G.l.c.(2) 

12 

99 

Minutes 

3-4 

FJRTIM(2,1) 

ft 

It 

ft 

ft 

5-6 

FIRTIM(3  ,1) 

*1 

ft 

ft 

tf 

7-8 

FIRTIM(4 ,1) 

tf 

It 

ft 

tt 

9-10 

FIRTIM (1,2) 

«f 

tf 

tf 

tt 

11-12 

FIRTIM (2, 2) 

If 

tf 

ft 

tt 

13-14 

FIRTIM(3,2) 

It 

ft 

ft 

tt 

15-16 

FIRTIM  (4, 2) 

tf 

ft 

tt 

ft 

17-18 

FIRTIM(1 , 3) 

ft 

ft 

ft 

ft 

19-20 

FIRTIM  (2, 3) 

ft 

ft 

tt 

ft 

21-22 

FIRTIM  (3, 3) 

tf 

tf 

ft 

ft 

23-24 

FIRTIM(4 ,3) 

ft 

ft 

ft 

ft 

25-26 

FIRTIM  (1,4) 

ft 

ft 

ft 

tf 

27-28 

FIRTIM (2, 4) 

«f 

ft 

tt 

tt 

29-30 

FIRTIM(3 ,4) 

If 

ft 

tt 

tt 

31-32 

FIRTIM (4, 4) 

ft 

ft 

t« 

tt 

Second 

Card:  Denominator  of 

Range  Component 

1-5 

DENOM(l ,1) 

G.l.c.(l) 

15 

99999 

Yards /Minute 

6-10 

DEN0M(2 ,1) 

tf 

i»*  ' 

tf 

ft 

(Same  order  of  subscripts  as  FIRTIM) 
•  • 

•  • 

76-80  DENOM(4,4)  "  "  " 

(IF  REPPOL  =  1,  skip  to  section  III.J. 7) 
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5 .  Refueling  Data 

a.  First  Card:  Miscellaneous 


Cols 

Vari able 

Reference 

Type 

Maximum 

Units 

1-4 

NOMO 

H.l.b 

14 

9999 

Minutes 

5-8 

COMBO 

H.3.C 

14 

9999 

Minutes 

(IF  REFSKD  =  2,  skip  to 

section  III 

.  J. .  5 .  c) 

b.  Next  4  x  FINFUL  Cards: 

User- 

Input  Refueling 

Schedule 

First 

Card 

1-7 

XFULSK (1,1,1) 

H.l 

F7.2 

9999.99 

Days 

8-14 

• 

• 

XFULSK(2 ,1,1) 

• 

»t 

ft 

ft 

ft 

• 

57-63 

• 

XFULSK (9, 1,1) 

ft 

tf 

tf 

tt 

Second 

Card 

1-7 

XFULSK(1 ,2,1) 

»t 

ft 

ft 

ft 

8-14 

• 

• 

XFULSK (2, 2,1) 

• 

•  t 

ft 

ft 

ft 

• 

57-63 

• 

XFULSK (9, 2,1) 

ft 

?» 

tt 

tf 

Third 

Card 

1-7 

XFULSK(1 ,3,1) 

ft 

ft 

tt 

tf 

8-14 

• 

• 

XFULSK (2, 3,1) 

■» 

ft 

Vf 

ft 

tt 

• 

57-63 

* 

XFULSK(9,3,1) 

ft 

ft 

ft 

♦f 

Fourth 

Card 

1-7 

XFULSK (1,4,1) 

ft 

ft 

ft 

tf 

8-14 

• 

• 

XFULSK(2',4,1) 

• 

• 

ft 

ft 

tf 

ft 

• 

57-63 

• 

XFULSK(9,4,1) 

ft 

ft 

ft 

tf 
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Next  4  Cards 


Same  as  first  four  cards  except  that  third  subscript  is  2. 


Continue  in  the  same  manner  until  there  are  FINFUL  sets  of 


4  cards 

each.  There 

must  be  4  x 

FINFUL 

cards  in 

the  user- 

input  refueling  schedule  section 

of  the 

data  deck. 

(Skip  to  section  III 

.J.S.d) 

c.  Next  Card:  Parameters  of  Computer- 

Calculated 

Refueling  Schedule 

Cols 

Variable 

Reference 

ZX££ 

Maximum 

Units 

1-6 

XCRFST(l) 

H.l.a 

F6.2 

999.99 

Days 

7-12 

XCRFST(2) 

ft 

tt 

Iff 

tt 

13-18 

XCRFST(3) 

t? 

tt 

»« 

ft 

19-24 

XCRFST(4) 

ft 

tt 

Iff 

tt 

25-30 

XCRFST(5) 

If 

tt 

*» 

tt 

31-36 

XDBFST(I) 

« 

91 

tt 

If 

fff 

37-42 

XDDFST(2) 

ft 

tt 

Iff 

I* 

43-48 

XDDFST(3) 

If 

ti 

If 

»« 

49-54 

XDDFST(4) 

tt 

ffff 

tt 

ft 

55-60 

XDDFST(5) 

Iff 

tt 

Iff 

tf 

61-65 

CRINT 

fl 

IS 

99999 

Minutes 

66-70 

DDINT 

It 

15 

99999 

Minutes 

d.  Next  Card:  Minimum 

and  Maximum  Refueling  Off 

Station  Times 

1*5 

MINFLC 

H.3.a 

IS 

99999 

Minutes 

6-10 

MAXFLC 

Iff 

tt 

ft 

If 

11-15 

MINFLD 

Iff 

tt 

tf 

II 

16-20 

MAXFLD 

Iff 

tt 

If 

tt 
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6.  Rearming  Data 

a.  first  Card:  Minimum  Allowed  Ammunition  Levels 
and  Minimum  and  Maximum  Rearming  Off-Station 
times  ”  & 


Cols 

Variable 

Reference 

1221 

Maximum 

Uni  > 

1-6 

MINLEV (1) 

H.  2 

16 

999999 

Rounds 

7-12 

MINLEV (2) 

It 

II 

II 

II 

13-18 

MINLF.V(3) 

It 

It 

II 

II 

19-24 

MINLEV (4) 

It 

•  1 

It 

If 

25-29 

M IN ARC 

H.3.b 

15 

99999 

Minutes 

30-34 

MAXARC 

ft 

It 

It 

II 

35-39 

MINARD 

It 

It 

II 

II 

40-44 

MAXARD 

It 

It 

It 

II 

b.  Next  4 

Cards : 

Initial  Ammunition 

Levels 

First 

Card 

1-7 

INAMM(1,1) 

H.2 

17 

9999999 

Rounds 

8-14 

• 

• 

INAMM(2 ,1) 

• 

It 

It 

ft 

II 

• 

57-63 

• 

INAMMC9 ,1) 

It 

It 

It 

II 

Second 

Card 

1-7 

INAMM(1 ,2) 

II 

It 

II 

II 

8-14 

INAMM(2 ,2) 

ft 

tt 

•  I 

It 

57-63  INAMM(9,2) 
Third  Card 


1-7  TNAMM(1 ,3) 

8-14  INAMM(2 ,3) 


II 


It 


57-63 


INAMM(9,3) 
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Fourth  Card 


1-7 

INAMM(1 ,4)  H . 2 

17 

9999999 

Rounds 

8-14 

INAMM(2,4)  M 

If 

It 

It 

57-63  INAMM(9 ,4) 

7.  Next  2  Cards:  Ammunition  Expenditures  Per  Minute 
Cols  Variable  Reference  Type  Maximum  Units 
First  Card 

1-4  MINAMM(l.l)  G.2  14  9999  Rounds/Minute 

5-8  MINAMM(2 ,1)  ••  "  "  " 

♦  • 

•  • 

(Same  order  of  subscripts  as  FIRTIM) 


61-64  MINAMM(4 ,4)  " 

Second  Card 

1-4  MAXAMM (1 ,1)  " 

5-8  MAXAMM(2 ,1)  " 


(Same  order  of  subscripts  as  FIRTIM) 


61-64  MAXAMM(4 ,4)  "  "  "  " 

8.  Last  Card:  Random  Number  Generator  Initializers 
This  card  contains  a  string  of  70  random  digits  in 
columns  1-70,  with  the  requirement  that  every  fifth  column 
contains  an  odd  number. 
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K.  COMPUTER  STORAGE  AND  TIME  REQUIREMENTS 

The  program  runs  on  an  IBM- 360/67  at  the  W.  R.  Church 
Computer  Center  at  the  U.  S.  Naval  Postgraduate  School.  For 
other  systems,  the  storage  and  time  requirements  may  differ 
from  those  stated  below. 

1 .  Computer  Storage 

The  amount  of  core  area  required  is  dependent  on  the 
amount  of  array  storage  set  aside  for  array  variables,  which 
of  course  depends  on  the  estimated  size  of  the  problem  to  be 
simulated.  The  program  as  presented  in  Appendix  B  allows  up 
to  2300  events  and  1000  missions  to  be  generated  during  each 
iteration  of  the  problem.  Thus  the  19  parameters  associated 
with  each  mission  all  have  dimension  size  1000  and  the  4  pa¬ 
rameters  associated  with  each  event  all  have  dimension  size 
2300.  With  these  array  sizes,  the  program  requires  approxi¬ 
mately  19SK  bytes  of  execution  core  area. 

2 .  Computer  Central  Processing  Unit  Time 

The  program  requires  approximately  1  minute  and  25 
seconds  of  CPU  time  for  the  compile  and  link  steps.  CPU  time 
required  in  execution  is  primarily  dependent  on  the  number 
of  missions  generated  and  the  number  of  iterations  performed. 
The  amount  of  execution  time  can  be  calculated  approximately 
as  1/200  seconds  per  mission  per  iteration.  The  expected 
number  of  missions  per  iteration  is  approximately^  XCFF(NR) . 

Therefore,  if  the  expected  number  of  missions  per 
iteration  for  a  problem  is  500  and  10  iterations  are  perform¬ 
ed  (NRIT  =  io),  then 

Execution  Time  =  500(1/200)  (10)  =  25  seconds. 
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IV.  VERIFICATION  AND  VALIDATION  OF  THE  MODEL 


A.  VERIFICATION 

The  fact  that  logic  flow  in  the  program  is  as  intended 
and  that  o*  ■'•out  statistics  accurately  describe  the  situa¬ 
tion  being  simulated  has  been  verified  in  several  ways. 

These  are  discussed  below. 

1 .  Verification  of  Intermediate  Calculations  and  Data 

Flow 

During  the  formulation  of  the  program  and  numerous 
debugging  runs,  results  of  calculations  and  other  pertinent 
data  at  various  points  in  the  program  were  printed  out  to 
determine  their  con*ectness .  When  the  program  was  considered 
finally  debugged,  such  checlcs  were  performed  at  all  points 
in  the  program  at  which  calculations  were  performed  or  sig¬ 
nificant  data  points  referenced.  The  results  of  these 
checks  were  carefully  examined,  and  all  values  were  deter¬ 
mined  to  be  correct. 

2 .  Tracing 

Logic  flow  through  all  possible  branches  of  the 
program  has  been  carefully  traced  and  determined  to  be  as 
intended. 

» 

3.  Extreme  Value  Problems 

When  the  program  was  considered  finally  debugged, 
several  problems  were  formulated  using  extreme  foreseeable 
maximum  and  minimum  values  of  input  parameters.  In  each 
case,  the  program  successfully  ran  to  completion  with 
statistical  results  reflecting  the  extreme  values. 
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.  Verification  of  Exponential  Distribution 

To  verify  that  mission  interarrival  times  are 
actually  produced  according  to  the  exponential  distribution, 
a  Kolmogorov-Smirnoff  test  for  goodness-of-fit  was  perform¬ 
ed  on  the  interarrival  times  in  a  subarea  picked  at  random 
in  a  trial  run  of  the  program.  The  subarea  picked  was  sub- 
area  11.  The  following  input  data  was  used  for  subarea  11 
in  the  trial  run: 

XCFF(ll)  =  35.0  missions 
XSHORT (11)  =  6.5  days 

XLONG(ll)  =  7.5  days 

Thus , 

XDURAT(ll)  =  7.0  days 

XCALRT(ll)  =  .00347  missions/minute. 

The  continuous  curve  in  Figure  3  is  a  plot  of  the 
CDF  of  an  exponential  distribution  with  rate  »  .00347.  The 
step  function  in  Figure  3  is  the  sample  CDF  of  the  23  mis¬ 
sion  interarrival  times  generated  for  subarea  11  in  the 
trial  run. 

The  Kolmogorov-Smirnoff  test  [pp.  238-240,  Ref.  5] 
consists  of  observing  the  absolute  difference  between  the 
value  of  the  sample  CDF  and  the  value  of  the  hypothesized 
CDF  at  all  points  in  the  sample.  The  test  statistic  is 

Dn  =  max  {  |  F(x)  -  Fn(x)  |  } 
n 

where  F(x)  is  the  value  of  the  hypothesized  CDF  at  x  and 
Fn(x)  is  the  observed  value  of  the  sample  CDF  at  x  for  a 
sample  of  size  n. 
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Figure  3\  Kolmogorov-Smirnoff  Test  -for  Goodness-of-Fit 

of  Mission  Interarrival  Times  to  the  Exponential 
Distribution. 


Inspection  of  Figure  3  yields  D„„  »  .165.  Table 

-  “  U  O 

VI  in  Reference  5  shows  that  D23  <  . 276  is  acceptable  with 
95%  confidence.  Therefore  the  sample  of  interarrival  times 
in  subarea  11  successfully  passes  the-  Kolmogorov- Smirnoff 
test  for  goodness-of-fit  for  the  exponential  distribution. 
From  this  result,  it  was  assumed  that  all  interarrival  times 
and  mission  firing  times  fit  an  exponential  distribution, 
since  the  same  method  of  generating  exponentially  distributed 
random  variables  was  used  in  each  case. 

5.  Comparison  with  GPSS  Model  - 

A  short  program  in  the  General  Purpose  Simulation 
System  (GPSS)  was  written  to  simulate  the  system  described 
herein  in  a  very  basic  mode,  which  was  easily  adaptable  to 
GPSS.  The  program  described  herein  was  run  using  input 
values  which  placed  it  in  exactly  the  same  mode  as  the  GPSS 
program.  The.  queueing  results  of  each  were  then  compared 
and  found  to  be  in  very  close  agreement. 

6.  Verification  of  Relationships  Between  Output 

Statistics 

By  investigating  the  mathematical  relationships  be¬ 
tween  various  output  statistics,  it  was  verified  that  the 
calculations  which  produce  the  output  values  are  properly 
performed.  An  example  of  an  investigation  of  such  a  rela¬ 
tionship  is  discussed  below. 

OVERALL  SHIP  UTILIZATION:  If  refueling  and  rearming 
are  not  used  in  the  problem,  overall  ship  utilization  should 
equal 
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* 


Total  Missions  x  Average  Service  Time  Per  Mission 
Number  of  Ships  x  Total  Problem  Time 

For  a  trial  run  of  the  program  in  which  refueling  and  rearm¬ 
ing  were  not  used,  the  following  results  were  obtained: 

NR855  =  NR647  =  NR538  =  1 
NR554  =  0 

TOTAL  NUMBER  OF  MISSIONS  FIRED  =  385 
AVERAGE  SERVICE  TIME  PER  MISSION  =43.0  minutes 
TOTAL  PROBLEM  TIME  =  8  days,  15  hours,  54  minutes 
OVERALL  SHIP  UTILIZATION  =  0.44. 

Using  these  figures  in  the  formula  stated  above  yields: 

Overall  Ship  _  _ 385  x  43.0 _  ... 

Utilization  3 x  (8  x  1440  +  IS  x  60  +  54) 

Since  overall  ship  utilization  is  not  calculated 
using  the  above  formula  in  the  program,  but  rather  by  accumu¬ 
lating  an  overall  sum  of  minutes  of  utilization  for  each 
ship,  this  result  lends  further  verification  to  both  logic 
flow  and  accuracy  of  all  of  the  above  listed  output  statis¬ 
tics.  By  methods  similar  to  the  above,  calculations  of 
other  related  output  statistics  were  also  verified. 

7.  Sensitivity  Testing 

As  a  final  method  of  verification,  the  sensitivity 
of  the  model  to  changes  in  input  parameters  was  tested.  An 
example  of  such  tests  is  discussed  below. 

WTDIS:  A  trial  run  of  the  program  was  made  using 
WTDIS  =  4000.  The  output  statistics  for  this  problem  yielded: 
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AVERAGE  TRAVEL  TIME  =16.5  minutes 

AVERAGE  WAITING  TIME  (ALL  MISSIONS)  =8.3  minutes 

PERCENT  OF  MISSIONS  REQUIRED  TO  WAIT  =  26.5%. 

The  same  problem  was  then  run  using  WTDIS  =  8000, 
all  other  input  values  remaining  unchanged.  It  was  ex¬ 
pected  that  each  of  these  output  statistics  would  increase 
by  amounts  estimated  as  follows: 

a.  Average  travel  time  should  increase  by  ap¬ 
proximately 

•70  |wr3T3j*  4-2  minutes 

since  each  ship  coming  from  the  waiting  pattern  would  have 
an  additional  4000  yards  to  travel  at  a  speed  of  20  knots, 
and  approximately  70%  of  the  missions  wore  fired  by  ships 
coming  from  the  waiting  pattern. 

b.  Average  waiting  time  (all  missions)  and  percent' 
of  missions  required  to  wait  should  increase  very  slightly 
due  to  the  small  increase  in  average  travel  time. 

The  results  of  the  trial  run  using  WTDIS  =  8000 
were  as  follows : 

AVERAGE  TRAVEL  TIME  =20.4  minutes 

AVERAGE  WAITING  TIME  (ALL  MISSIONS)  =8.5  minutes 

PERCENT  OF  MISSIONS  REQUIRED  TO  WAIT  =27.2 

These  results  are  very  close  to  those  expected, 
further  verifying  logic  flow  and  program  results. 
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B.  VALIDATION 

The  simulation  cannot  be  considered  validated  since  it 
is  obviously  infeasible  to  actually  physically  conduct  a 
problem  for  experimental  purposes  within  the  scope  of  this 
thesis,  and  no  experimental  data  of  similar  exercises  is 
available.  The  design  of  the  model  is  primarily  based  on 
personal  experience  in  NGFS  both  from  aboard  ship  and  as  a 
spotter  ashore,  and  the  results  obtained  by  the  simulation 
model  appear  to  be  valid  within  the  limits  of  that  experience. 
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V.  PARTIAL  STDF  FEASIBILITY  STUDY 


j 

! 


With  the  computer  simulation  model  finalized  and  fully 
tested,  a  partial  study  of  the  feasibility  of  the  STDF  con¬ 
cept  was  undertaken.  This  study  focused  only  upon  the  que¬ 
ueing  aspects  of  providing  NGFS  from  the  STDF  to  SIAF's 
ashore. 

A.  OBJECTIVE  OF  THE  STUDY 

The  overall  objective  of  the  study  can  be  generally 
stated  as  follows: 

’’Determine  the  optimum  number  and  mix  of  NGFS  ships 
required  to  adequately  support  SIAF’s  in  an  insur¬ 
gent-contested  area.” 

Before  proceeding  further,  the  term  "adequate”  support 
must  be  defined.  Since  the  model  does  not  consider  results 
of  missions  fired,  adequate  support  cannot  be  defined  in 
terms  of  results  achieved.  Rather,  the  definition  must  be 
in  terms -of  the  queueing  aspects  of  NGFS. 

1.  Decision  Criteria 

Two  arbitrary  decision  criteria  were  established 
in  order  to  define  ’’adequate”  support.  These  criteria  are 
based  solely  on  personal  experience  and  cannot  be  considered 
authoritative.  The  decision  criteria  are: 

a.  The  average  waiting  time  per  mission  should  be 
less  than  10  minutes.  Spotting  teams  ashore  become  quite 
anxious  and  tneir  position  becomes  increasingly  tenuous  as 
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the  length  of  time  between  their  call  for  fiTe  and  notifica¬ 
tion  that  a  ship  is  on  the  way  increases.  Ten  minutes  seems 
like  a  reasonable  upper  bound  on  the  length  of  time  that  a 
patrol  should  have  to  expect  to  wait. 

b.  Overall  ship  utilization  should  be  at  least 
SOI.  Paucity  of  ships  to  assign  to  a  particular  operation 
dictates  a  requirement  such  as  this ,  and  is  the  basis  of  the 
need  for  a  feasibility  study  of  the  queueing  aspects  of  NGFS 
of  SIAF*s . 

2.  Objective  Function 

Having  established  the  two  decision  criteria  to 
define  "adequate**  support,  the  objective  of  the  study  be¬ 
comes  : 

"Determine  the  minimum  number  (and  mix)  of  NGFS 
ships  required  to  satisfy  the  decision  criteria.’* 

In  terms  of  the  variables  and  the  output  statistics 
in  the  computer  model,  the  objective  can  be  stated  as: 

Minimize:  NR8SS  +  NR647  ♦  NR554  +  NR538 
Subject  to:  AVERAGE  WAITING  TIME  (ALL  MISSIONS)  <10.0 
OVERALL  SHIP  UTILIZATION  >  0.S0. 

3.  Optimization  Technique 

The  optimization  technique  used  can  be  described 

generally  as  a  3-step  procedure.  These  steps  are: 

a.  Hypothesize  a  typical  STDF/SIAF  situation. 

b.  Exercise  the  simulation  model,  varying  the 


number  of  ships  of  each  type  which  are  present.- 


c.  Inspect  the  results  of  the  simulation  exercise 
to  determine  the  minimum  number  (and  mix)  of  ships  which 
best  satisfies  the  decision  criteria. 

B.  HYPOTHESIZED  SITUATION  -  SCENARIO 

In  order  to  exercise  the  computer  simulation,  a  scenario 
was  developed  which  seems  to  describe  a  typical  situation  as 
envisioned  under  the  STDF/SIAF  concept. 

1 .  General  Situation 

A  small  country,  allied  with  the  U.S.,  is  under 
imminent  threat  of  insurgent  takeover.  For  the  past  several 
months,  there  have  been  numerous  reports  of  infiltration  into 
the  country  and  tensions  have  increased  with  the  country's 
unfriendly  neighbor  to  the  southwest.  During  the  past  two 
weeks  there  have  been  terrorist  attacks  in  the  small  fishing 
villages  located  along  the  banks  of  a  river,  which  divides 
the  country  roughly  in  half.  The  capitol  city  is  located  at 
the  mouth  of  this  river.  The  government  now  believes  that 
the  infiltrators  are  an  insurgent  force  from  the  unfriendly 
neighboring  country,  and  that  the  insurgents  are  building  up 
for  an  attack  on  the  capitol  city  in  an  attempt  to  unseat 
the  present  government.  The  allie^  .country  has  requested 
emergency  military  assistance  from  the  U.S.  to  suppress  the 
insurgent  force. 

2 .  Enemy  Situation 

The  recent  terrorist  attacks  and  other  reports  in- 
dicate  that  the  bulk  of  the  enemy  force  is  dispersed  along 
the  banks  of  the  river  in  small  guerilla  groups.  However, 
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sines  there  have  been  no  sightings  or  use  of  heavy  artillery 
or  other  significantly  large  military  equipment  near  the 
river,  it  is  believed  that  some  of  the  enemy  force  is  en¬ 
gaged  in  transporting  heavy  equipment  to  strategic  locations 
to  the  west  of  the  capitol  city.  The  most  likely  routes  for 
movement  of  such  equipment  are  through  the  sparsely  populated 
mountainous  regions  in  the  southwest  and  northwest  parts  of 
the  country.  There  are  estimated  to  be  between  7500  and 
10000  insurgents  within  the  country  at  the  present  time. 

3.  Friendly  Situation 

The  U.S.  has  agreed  to  send  elements  of  the  STDF 
to  assist  in  suppressing  the  insurgent  force.  However,  the 
U.S.  has  made  it  clear  that  this  military  assistance  is  on 
an  emergency  basis  only,  and  will  be  provided  for  a  period 
of  seven  days  only.  The  STDF  commander  has  received  specific 
orders  to  place  only  a  minimum  number  of  men  ashore,  and  that 
patrols  placed  ashore  are  not  to  actually  engage  in  combat 
except  in  self-defense. 

a.  NGFS/SIAF  Detachment 

The  STDF  commander  detached  the  following  ships 
from  the  STDF  to  provide  the  emergency  assistance  requested: 

1  LHA  (with  one  Marine  Expeditionary  Unit 
embarked) 

1  AOE 

X  Cruisers  and  Destroyers 

b.  Area  of  Operations 

(1)  Objective  Area.  The  objective  area  is  de¬ 
fined  as  the  area  surrounding  the  capitol  city,  10  miles  north 
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and  south  of  the  river,  and  including  up  to  the  peaks  of  the 
mountain  ranges  in  the  northwest  and  southwest.  The  objective 
area  is  approximately  20  miles  wide  and  10  miles  deep  (inland 
from  the  coast) . 

(2)  Patrol  Subareas.  Based  on  expected  densi¬ 
ties  of  enemy  targets  and  terrain  features,  the  objective  area 
has  been  divided  into  34  SIAF  patrol  subareas.  Patrols  will 
be  inserted  into  these  subareas  by  helicopter  beginning  at 
0500  (H-hour)  on  the  first  day.  Order  and  time  of  patrol  in¬ 
sertion  will  be  randomly  picked.  All  patrols  will  have  been 
inserted  no  later  than  2100  on  the  first  day. 

(3)  Waiting  Pattern.  When  not  assigned  to  fire 
a  mission,  ships  will  operate  within  a  10  mile  wide  area  4-6 
miles  offshore  and  centered  on  the  center  of  the  objective 
area  since  the  heaviest  target  densities  are  expected  to  be 
near  the  center  of  the  objective  area. 

c.  Firing  Doctrine 

In  this  situation,  the  U.S.  is  very  much  con¬ 
cerned  that  its  troops  ashore  not  engage  in  combat  unless 
absolutely  necessary.  Thus  it  is  desired  that  all  calls  for 
fire  from  SIAF  patrols  be  answered  as  rapidly  as  possible  by 
any  ship  available,  and  that  personnel  targets  be  fired  at 
before  material  targets  if  both  occur  simultaneously. 

d.  Refueling  and  Rearming 

(1)  Location  of  APE.  The  assigned  AOE  will  re¬ 
main  approximately  15  miles  offshore  to  accomplish  refueling 
and  rearming.  Thus,  a  one-way  trip  from  the  waiting  pattern 
to  the  AOE  at  a  speed  of  20  knots  should  require  approximately 
30  minutes. 
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(2)  Refueling  Schedules.  Since  it  is  antici¬ 


pated  that  the  operation  will  last  approximately  one  week,  no 
ship  should  require  refueling  more  than  once.  Destroyers  are 
scheduled  to  refuel  beginning  at  1100  on  the  fifth  day,  with 
other  destroyers  following  the  first  at  four-hour  intervals. 
Destroyers  may  also  refuel  beginning  at  1100  on  the  ninth  day 
if  the  operation  lasts  that  long.  Cruisers  are  scheduled  to 
refuel  beginning  at  1100  on  the  sixth  day,  with  other  cruisers 
following  the  first  at  four-hour  intervals.  It  is  anticipated 
that  destroyers  should  require  60  to  90  minutes  alongside  the 
AOE  during  each  refueling,  while  cruisers  should  require  90 
to  120  minutes.  As  the  operation  progresses  and  the  time  of 
withdrawal  of  all  patrols  from  the  objective  area  becomes 
reasonably  firm,  no  ship  will  depart  station  for  refueling 
within  12  hours  of  this  anticipated  time  of  total  withdrawal. 

(3)  Initial  and  Minimum  Ammo  Levels.  Prior  to 
arrival  at  the  area  of  operations,  all  NGFS  ships  have  topped 
off  their  onboard  ammunition  levels.  The  S7DF  commander  has 
prescribed  minimum  allowable  ammunition  levels  for  each  type 
of  ship.  These  ammunition  levels  are  as  follows: 

Type  Ship  Initial  Minimum 


8"/55 
6"/47 
5’754 
5'*/  38 


8000  .  . 

1000 

10000  1 

1000 

12000 

2000 

6000 

1000 

It  is  anticipated  that  destroyers  should  require  2  to  3  hours 
alongside  the  AOE  to  rearm,  while  cruisers  should  require  3 
to  4  hours. 
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C.  HYPOTHESIZED  SITUATION  -  INPUT  VARIABLES 


Having  hypothesized  what  appears  to  be  a  typical  STDF/ 
SIAF  situation,  the  next  step  was  to  translate  the  scenario 
into  values  of  input  variables  in  order  to  exercise  the  com¬ 
puter  simulation.  Appendix  C  is  the  printout  of  input  data 
used  in  the  hypothesized  problem,  for  a  run  in  which  NR85S  = 
NR647  =  NR554  =  1  and  NR538  =  0. 

1.  Zero  Problem  lime 

Zero  problem  time  is  0500  on  the  first  day  of  the 
operation.  Thus,  HHOUR  =  300.  All  insertion  times  are  in 
minutes  after  0500,  i.e.,  INSERT(2)  ■  140  indicates  a  SIAF 
patrol  will  be  -nserted  into  subarea  2  at  0720  on  the  first 
day. 

2.  Problem  Area  Map 

Figure  4  is  a  map  of  the  problem  area  translated  into 
the  coordinate  system  described  in  Section  III.B.  Subarea 
boundaries  used  are  not  necessarily  exactly  as  they  would  be 
in  the  actual  operation,  i.e,,  the  river  would  probably  be 
the  boundary  of  patrol  subareas  on  either  side  of  it.  The 
following  input  information  is  shown  on  the  map  for  each  sub- 
area: 


NR 

XCFF(NR) 

INSERT  (NR) 

XPR0B1  (NR) 
XPR0B2 (NR) 
XP  ROB  3  (NR) 
XPR0B4 (NR) 

The  values  used  for  expected  number  of  calls  for 
fire  and  target  type  probabilities  are  based  on  the  expected 
locations,  density,  and  type  of  enemy  targets  as  indicated 
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Figure  4.  Problem  Area  Map  for  STDF  Feasibility  Study. 


in  Section  V.B. 2.  XCFF  and  XPR0B1  +  XPR0B2  are  high  near  the 
river,  and  XPR0B2  >  XPR0B1  since  it  is  believed  that  the  in¬ 
surgent  force  operates  primarily  in  small  guerilla  groups. 
XPR0B3  +  XPR0B4  is  high  in  the  mountainous  areas  to  the 
southwest  and  northwest,  indicating  suspected  movement  of 
supplies  through  those  areas,  with  XPR0B3  >  XPR0B4  indicat¬ 
ing  primarily  heavy  equipment  being  moved. 

3.  Expected  Patrol  Durations 

The  maximum  and  minimum  patrol  durations  were 
chosen  such  that  the  expected  time  of  extraction  of  each 
patrol  should  be  approximately  7  days  after  zero  problem 
time.  However,  patrol  duration  parameters  are  given  a  wider 
range  of  values  in  the  areas  of  heavy  density  of  enemy  per¬ 
sonnel.  This  reflects  a  higher  probability  of  early  extrac¬ 
tion  due  to  SIAF  personnel  casualties. 

4.  Assignment  Doctrine  and  Queue  Discipline 

The  discussion  in  Section ‘V.B. 3. c  implies  that  ‘’any 
ship  with  guns  of  adequate  range  may  fire"  any  mission  (ASNDOC 
-  1) ,  and  "personnel  missions  have  priority"  in  the  queue 
(QDISC  =  2). 

5.  Refueling  Schedule 

The  refueling  schedule  as 'discussed  in  Section  V.B. 

3. a  translates  into  the  following  values  of  refueling  sched¬ 
ule  variables: 

REPPOL  *  2 

REFSKD  =  1  (Computer  calculated  schedule  used) 

COMBO  =  60  (Time  for  round  trip  between  AOE  and 
waiting  pattern  is  one  hour) 
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NOMO  =  720  (No  more  refuelings  within  12  hours  of 
latest  extraction  time) 

XCRFST(l)  ■  5.2S  (S  days,  6  hours  after  zero  prob¬ 
lem  time  is  1100  on  6th  day) 

XCRFST(2)  thru  XCRFST(S)  =  99.00 

XDDFST(l)  =4.25 

XDDFST(2)  =8.25 

XDDFST(3)  thru  XDDFST(5)  =  99.00 

CRINT  =  240 

DDINT  =  240 

Other  values  of  refueling  and  rearming  parameters 
are  translated  as  indicated  in  Appendix  C. 

6.  Firing  Time  Parameters 

The  components  of  mean  firing  time  as  indicated  in 
Appendix  C  were  determined  as  follows: 

a.  Denominator  of  Range  Factor 

The  range  to  the  target  generally  has  less  ef¬ 
fect  on  the  accuracy  of  8"/S5  and  6"/47  shells  than  5"  shells 
since  they  are  larger  and  designed  for  longer  ranges. 

b.  Mean  Firing  Time  Due  to  Target  Type 

A  5"/54  gun  is  considerably  more  effective 
against  personnel  targets  and  should  require  considerably 
less  time  to  fire  a  personnel  target  mission  due  to  its 
rapid  rate  of  fire.  8"/55  and  6”/47  guns  require  less  time 
against  material  targets  due  to  their  high  explosive  power, 
while  a  5”/54  gun  requires  less  time  than  a  5"/38  due  to  its 
rapid  rate  of  fire. 
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D.  INVESTIGATION  OF  VARIOUS  NUMBERS  AND  MIXES  OF  SHTPS 

PRESENT 

1.  Two  Ships 

Rough  calculations  based  on  expected  number  of  mis¬ 
sions,  expected  service  time,  and  expected  problem  duration 
indicate  that  it  is  very  unlikely  that  any  combination  of 
two  ships  would  be  able  to  provide  adequate  support.  These 
rough  calculations  are  as  follows: 

34 

Y,  XCFF(NR)  =  414.0 
NR“1 

Expected  total  service  time  per  mission  s  45  minutes 

Expected  problem  duration  s  7  days 

414  missions  x  3/4  hours/mission  s  310  ship  hours 
required 

Two  ships  x  7  days  x  24  hours/day  *  336  ship  hours 
available. 

The  variability  of  mission  interarrival  times,  mis¬ 
sion  firing  times,  and  travel  distances,  and  the  requirements 
for  refueling  and  rearming  would  most  likely  cause  an  Unac¬ 
ceptable  queueing  situation  to  develop.  Therefore,  combina¬ 
tions  of  two  ships  were  not  investigated  initially. 

2.  Three  Ships 

It  is  not  at  all  apparent  from  rough  calculations 
whether  or  not  three  ships  would  be  able  to  adequately  sup¬ 
port  the  hypothesized  operation.  Herein  lies  the  basic 
justification  for  the  use  of  simulation,  since  simulation 
accounts  for  the  variability  in  problem  parameters  and  pro¬ 
duces  an  average  (expected)  result  on  which  to  base  decisions. 
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To  investigate  the  queueing  results  for  three  ships 
present,  the  problem  was  run  16  times,  representing  all  pos¬ 
sible  combinations  of  0  to  three  ships  of  each  type  which 
add  to  3,  with  the  requirement  that  at  least  one  8"/55  ship 
or  one  5"/54  ship  be  present,  since  the  northern  extremity 
of  the  objective  area  is  25000  yards  from  the  gunline.  Ten 
iterations  per  run  (NRIT  =  10)  were  used  to  determine  aver¬ 
aged  results  and  standard  deviations.  Table  I,  on  the  next 
page,  presents  the  two  decision  criteria  output  statistics 
for  each  of  the  16  combinations  of  ships.  Figures  in  the 
table  are  the  upper  and  lower  limits  of  a  951  confidence 
interval.  These  limits  were  calculated  using  the  equation 

x  -  2.26  0=  <  n  <  x  +  2.26  o* 

which  was  derived  as  discussed  in  Section  III. I. 4. 

Inspection  of  the  figures  in  Table  I  shows  that  re¬ 
gardless  of  the  mix  of  ships,  overall  ship  utilization  is 
approximately  50%  when  three  ships  are  present.  From  this 
observation  it  can  be  readily  deduced  that  utilization  when 
four  ships  are  present  would  be  considerably  less  (probably 
about  35%)  while  utilization  when  only  two  ships  are  present 
would  be  greater  (probably  about  63%) .  Based  on  this  con¬ 
clusion,  the  simulation  was  exercised  for  the  seven  possible 
combinations  of  two  ships  to  see  if  any  combination  would 
yield  a  satisfactory  average  waiting  time.  The  results  of 
these  seven  runs  are  presented  in  Table  II,  on  the  next 
page. 
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TABLE  I.  Decision  Criteria  Output  Statistics  for  Combinations 


0 

0 

0 

0 

0 

0 

1 

1 

1 

X 

1 

1 

2 

2 

2 

3 


NR647 

NR554 

NR538 

AVG  WAITING  TIME 
MIN  MAX 

wwtwAwu  utur 

UTILIZATION 

MIN  MAX 

0 

1 

2 

24.0 

35.0 

.482 

.528  j 

0 

2 

1 

.  7.0 

10.2 

.484 

.510  1 

0 

3 

0 

3.7 

5.3 

.447 

.485 

1 

1 

1 

8.6 

13.4 

.470 

.522 

1 

2 

0 

4.0 

6.2 

.434 

.494  j 

2 

1 

0 

5.6 

7.8 

.477 

.493  | 

0 

0 

2 

30.5 

38.1 

.487 

•  539  j 

0 

1 

1 

7.6 

10.2 

.476 

.516  j 

0 

2 

0 

3.6 

6.6 

.460 

.508  .  j 

1 

0 

1 

7.8 

11.2 

.478 

.524  j 

1 

1 

0 

3.7 

5.7 

.459 

.481  1 

2 

0 

0 

1 

0 


0 

0 

1 

0 

0 


0 

1 

0 

0 

0 


6.4 

8.5 

4.6 

4.7 
4.5 


8.0 

11.1 

8.4 

7.3 

9.9 


.475 

.494 

.478 

.476 

.468 


.507 

.526 

.518 

.522 

.508 


TABLE  II. 


ofCTwo°ShiJsteria  °UtpUt  Statistics  for  Combinations 


NR855  NR647  NR554  NR538 


0 

0 

0 

1 

1 

1 

2 


AVG  WAITING  TIME 
MIN  MAX 


0 

0 

1 

0 

0 

1 

0 


1 

2 

1 

0 

1 

0 

0 


OVERALL  SHIP 
UTILIZATION 
MIN  MAX 


1 

0 

0 

1 

0 

0 

0 


48.8 

73.4 

.687 

.743 

26.0 

41.6 

.676 

.742 

29.8 

42.8 

.679 

.715 

65.5 

85.9 

.719 

.767 

25.6 

44.4 

.666 

.732 

33.4 

45.8 

.674 

.743 

33.1 

46.9 

.685 

.751 
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At  this  point  in  the  analysis,  it  was  observed  that 
no  combination  of  ships  present  would  simultaneously  satisfy 
the  two  decision  criteria  within  the  limits  of  a  951  con¬ 
fidence  interval.  When  three  ships  are  present,  there  is  no 
combination  for  which  it  can  be  said  with  95%  confidence  that 
OVERALL  SHIP  UTILIZATION  >  .50,  while  there  is  no  combination 
of  two  ships  for  which  it  can  be  said  with  95%  confidence 
that  AVERAGE  WAITING  TIME  <  10.0  minutes.  Thus,  one  of  the 
two  decision  criteria  had  to  be  relaxed.  Inspection  of 
Tables  I  and  II  made  it  readily  obvious  that  the  ship  uti¬ 
lization  criteria  was  the  logical  criteria  to  relax.  Thus, 
the  decision  criteria  was  changed  to  OVERALL  SHIP  UTILIZA¬ 
TION  >  .43. 


With  these  revised  decision  criteria,  the  following 
nine  combinations  of  three  ships  satisfy  both  criteria  simul¬ 
taneously: 


NR855 

0 

.  0 
0 
1 
1 
1 
2 
2 
3 


NR647 

0 

1 

2 

0 

1 

2 

0 

1 

0 


NR5S4 

3 

2 

1 

2 

1 

0 

1 

0 

0 


NR538 

0 

0 

0 

0 

0 

0 

0 

0 

0 


The  striking  characteristic  of  these  nine  combina¬ 
tions  is  that  NR538  *  0  in  all  of  them  and  they  are  the  only 


89 


nil*,  combinations  in  which  NR538  *  0.  This  is  obviously  be¬ 
cause  of  the  short  range  of  the  5"/38  gun,  which  limits  its 
usability  when  the  objective  area  is  more  than  7  miles  deep. 

In  an  attempt  to  further  distinguish  between  these 
nine  combinations  of  three  ships,  two  other  output  statistics 
were  investigated.  These  were  MAXIMUM  WAITING  TIME  and  PER¬ 
CENT  OF  MISSIONS  REQUIRED  TO  WAIT.  Table  III  on  the  next 
page  presents  these  statistics  for  the  nine  combinations  of 
ships  under  consideration. 

Inspection  of  Table  III  indicates  that  two  combina¬ 
tions  are  considerably  less  desirable  than  the  others. 

These  are: 


NR8S5 

NR647 

NR554 

0 

2 

1 

1 

2 

0 

"he  striking  characteristic  of  these  combinations 
is  that  NR647  =  2  in  each  of  them,  and  they  are  the  only  two 
combinations  of  the  original  16  in  which  NR647  =  2. 

E.  CONCLUSION  OF  THE  STUDY 

Based  on  the  foregoing  analysis  of  the  output  statistics, 
the  fol/'-wing  conclusion  was  reached: 

"For  u  .ituaticn  such  as  described  herein,  any  combina¬ 
tion  of  three  NGFS  ships  of  which  none  are  5"/38  ships  and 
less  than  two  are  6"/47  ships  is  optimal  in  terms  of  providing 
adequate  support  to  SIAF’s  ashore  in  an  insurgent-contested 
area." 
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Table  III.  MAXIMUM  WAITING  TIME  and  PERCENT  REQUIRED  TO  WAIT 
Statistics  for  the  Nine  Best  Combinations  of 
Three  Ships 

• 

PERCENT 

MAX  WAITING  TIME  REQUIRED  TO  WAIT 


NR855 

NR647 

NR554 

MIN 

MAX 

MIN 

MAX 

0 

0 

3 

69 

141 

15.5 

32.2 

0 

1 

2 

75 

169 

13.6 

36.6 

0 

2 

1 

183 

251 

21.8 

31.8 

1 

0 

2 

55 

123 

13.7 

39.7 

1 

1 

1 

71 

111 

13.8 

34.8 

1 

2 

0 

184 

342 

21.6 

33.8 

2 

0 

1 

80 

220 

18.3 

37.7 

2 

1 

0 

88 

174 

18.2 

37.2 

3 

0 

0 

110 

254 

11.7 

41.7 

I 

1 

=1 

-i 
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In  terms  of  STDF  feasibility,  this  can  be  interpreted 
to  mean; 

"If  three  NGFS  ships,  of  which  none  are  5"/38  ships  and 
less  than  two  are  6"/47  ships,  are  available  for  inclusion 
in  the  STDF,  then  the  STDF  concept  is  feasible  in  terms  of 
providing  NGFS  to  SIAF’s  ashore  in  a  situation  such  as 
described  herein. " 

It  must  be  emphasized  that  these  conclusions  only  apply 
to  the  hypothesized  situation.  The  sensitivity  of  the  con¬ 
clusions  to  changes  in  input  parameters  has  not  been  tested 
as  would  certainly  be  necessary  in  a  detailed  feasibility 
study.  It  is  not  considered  within  the  scope  of  this  thesis 
to  undertake  such  an  extensive  study.  Rather,  it  is  con¬ 
sidered  that  the  objectives  stated  in  Section  I.C  have  been 
met  at  this  point. 
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VI;  CONCLUSIONS 

It  is  recognized  that  there  are  several  ways  in  which 
the  model  presented  in  this  paper  is  capable  of  refinement 
and  extension  in  order  to  more  accurately  simulate  realistic 
SIAF  support  situations  and  be  more  useful  in  investigating 
the  overall  feasibility  of  the  STDF  concept.  Although  such 
refinement  and  extensions  were  not  considered  to  be  within 
the  scope  of  this  thesis,  some  of  them  are  pointed  out  in 
this  section  for  possible  future  consideration  by  users  of 
the  model. 

A.  REFINEMENT  OF  THE  MODEL 

1 .  Travel  Distance  Calculations 

The  calculations  of  travel  distances  on  which  ship 
assignments  are  based  can  be  made  more  realistic  by  develop¬ 
ing  a  dynamic  scheme  of  keeping  track  of  actual  positions  of 
unassigned  ships  throughout  the  run  of  the  model. 

2 .  Ship  Assignment  Doctrines 

Ship  assignment  doctrines  other  than  the  two  used 
in  this  model  can  be  devised  to  simulate  a  wider  range  of 

realistic  options  available  to  the  STDF  commander  in  as- 

*  • 

signing  ships  to  missions. 

3.  Queue  Discipline 

A  preempt  option  for  ’’urgent”  missions  can  be  in¬ 
troduced  to  add  realism  to  the  model. 

4.  Refueling  and  Rearming  Off-Station  Times 

The  amount  of  time  spent  by  each  ship  off-station 
for  rearming  and/or  refueling  can  be  made  more  realistic  by 
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developing  a  method  of  simulating  movement  of  oilers  and  am¬ 
munition  ships  throughout  the  problem. 

5 .  SIAF  Patrol  Movements 

The  model  can  be  made  more  realistic  by  developing 
a  dynamic  scheme  of  simulating  movement  of  SIAF  patrols  with¬ 
in  their  assigned  subareas.  Such  a  scheme  could  possibly 
be  the  linking  of  this  model  with  the  TRW  SIAF  model  in  some 
way. 

6 .  Distributions  of  Variable  Quantities 

All  of  the  distributions  of  variable  quantities 
used  in  the  model  are  based  primarily  on  personal  experience 
in  Naval  operations.  It  is  possible  that  some  of  the  dis¬ 
tributions  do  not  reflect  realistic  situations  as  accurately 
as  some  other  distribution  might.  Refinement  of  the  model 
should  include  an  investigation  of  historical  or  experimental 
data  to  determine  the  actual  distributions  of  the  variable 
quantities.  For  example,  data  concerning  actual  mission 
firing  times  should  be  analyzed  to  determine  if  firing  times 
actually  conform  to  the  exponential  distribution  or  to  some 
other  probability  distribution. 

B.  EXTENSIONS  OF  THE  MODEL 

1.  Inclusion  of  Air  Strikes 

The  model  can  be  extended  to  simulate  air  strikes  by 
helicopters  or  fixed-wing  aircraft  originating  from  carriers 
assigned  to  the  STDF  when  targets  could  be  more  effectively 
encountered  by  air  strikes  or  are  beyond  the  maximum  range  of 
Naval  gunfire. 
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2 .  Inclusion  of  Battleships 

Although  there  are  no  battleships  in  commission  at 
the  present  time,  the  possibility  that  battleships  could  be 
brought  back  into  use  in  the  future,  as  was  the  case  of  the 
NEW  JERSEY  in  the  Vietnam  War,  indicates  the  desirability  of 
including  battleships  in  the  model.  However,  if  no  ships  of 
a  particular  type  were  present  in  a  particular  problem,  the 
model  in  its  present  state  could  be  used  to  simulate  the 
presence  of  battleships  by  assigning  battleship  parameters 
to  the  input  variables  normally  used  for  ships  of  the  type 
not  present.  This  might  require  some  modification  to  the 
ship  assignment  doctrine. 

3.  Multiple  Patrols  in  a  Subarea 

The  model  could  be  extended  to  include  the  option 
of  inserting  a  second  patrol  into  a  subarea  when  the  first 
patrol  is  extracted,  and  so  on. 

C.  OVERALL  CONCLUSION 

While  it  is  recognized  that  the  model  presented  in  this 
paper  is  capable  of  refinement  and  extension  in  several  ways, 
it  is  considered  that  the  model  in  its  present  state  is  use¬ 
ful  as  a  means  of  providing  general^  ideas  concerning  the 
queueing  aspects  of  a  wide  variety  of  user-defined  NGFS/SIAF 
situations.  It  is  also  considered  that  the  model  can  be 
useful  in  operational  planning,  i.e.,  the  model  could  be  used 
to  investigate  the  effects  of  locating  the  waiting  pattern  in 
various  areas,  etc. 
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COMPUTER  PROGRAM  GENERALIZED  FLOW  CHART 


r 


I 


CALCULATE  TARGET  COORDINATES  AND  TARGET  TYPE  j 


DETERMINE  MISSION 
TYPE  BASED  ON  TARGET 
RANGE  ONLY 


DETERMINE  MISSION 
TYPE  BASED  ON  RANGE 
AND  TARGET  TYPE 


CHECK  AVAILABILITY  OF  SHIPS  ELIGIBLE 
TO  FIRE  THIS  TYPE  OF  MISSION 


CALCULATE  TRAVEL 
DISTANCE  OF  EACH 
AVAILABLE  SHIP, 
ASSIGN  CLOSEST  SHIP 
70  FIRE  THE  MISSION, 
SET  ASSIGNED  SHIP 
IN  A  FIRING  STATUS 


I 


CALCULATE  SERVICE 
TIME  AND  MISSION 
COMPLETION  TIME, 
AND  PLACE  MISSION 
COMPLETION  ON  THE 
EVENT  CHAIN 


PUT  MISSION 
IS  THE  QUEUE 
AHEAD  OF  ALL 
OTHERS  WITH 
TARGET  TYPES 
3  OR  4 


PUT  MISSION 
AT  END  OF 
THE  QUEUE 


1=1  +  1 


0 
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98 
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ASSIGN  FIRST  ELIGIBLE  MISSION, 
I?  ANY,  IN  THE  QUEUE  TO  BE 
FIRED  BY  THIS  SHIP.  SET 
SHI?  IN  A  "FIRING"  STATUS 


t 

CALCULATE  SERVICE  TIKE  AND 
MISSION  COMPLETION  TIME  OF 
ASSIGNED,  AND  PLACE  MISSION 
COMPLETION  ON  THE  EVENT  CHAIN 

* 

MOVE  ALL  SUCCEEDING  MISSIONS 
UF  ONE  POSITION  IN  THE  QUEUE 


APPENDIX  B:  COMPLETE  COMPUTER  PROGRAM  LISTING 
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APPENDIX  C:  PRINTOUT  OF  INPUT  DATA  FOR  STDF  FEASIBILITY  STUDY 
(For  the  Case  When  NR85S  =  NRSS4  *  1,  NR533  »  0) 


***  INPUT  RATA  *** 


NUMBER  0p  8*»/55  SHIPS  °RESCNT  *  1 
MIS  SI  SHIPS  PRESENT  »  1 
NUMBER  OF  5"/5 4  SHIPS  PPESj-NT  =  I 
NUMBER  OP  5"/38  SHIPS  PREScNT  =  0 


MAX  RANGE  OF  8V.35  G;  N  *  30000  VOS 
MAX  RANGE  OF  6**/47  Go. I  *  21000  VOS 
MAX  RANGE  OF  S-/54  GUN  =  26000  YDS 
MAX  RANGE  OF  5«/38  GUN  *  15000  YDS 


nrl-riv?!  SI  ££§J£52.,ES£le°l,.MAinNG  PATTERN  FROM  Y-AXIS  »  10000  YDS 

i^LENiTHF0P05AlflNGCPAfTi=N  »I20000PYDS€RN  FR°"  6UNtI*5  *  5000  Y°S 

N-S  LENGTH  OF  WAITING  PATTERN  *  5000  YOS 

S?sE§6?8pFttJ1G.PO!*3ilS2  * 20  KN0TS 

MAX  SETUP  TI*»E  *  15  MINS 

ASSIGNMENT  DOCTRINE:  1 

0UEU6  DISCIPLINE:  2 

REPLENISHMENT  POLICY:  2 

TYPE  OF  FUELING  SCHEDULE  UScO:  2 


*  COMPONENTS  OF  MEAN  FIRING  TIME  * 

•*EAN  FIRING  TIME  DUE  TO  TARGET  TYPE:  (MINUTES  FRCM  TABLE  BELOW) 


TYPE  SHIP 


TARCET  TYPE 
12  3  4 


15 

10 

12 

10 

14 

8 

11 

9 

8 

5 

14 

10 

18 

12 

18 

12 

DENOMINATOR  OF  RANGE  FACTOR:  (FROM  TABLE  BELOW) 


TYPE  SHIP 


TARGET  TYPE 
2  3 


10000  8000  6000  6000 

100CC  8000  6000  6000 

10000  8000  6000  6000 

8000  4000  4000  4000 
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*nc»  inn.'*”&T  ICN  * 


25000 

2500C 

25000 

250CC 

25000 

25000 

2500C 

25000 

20000 

200C0 

20000 

20000 

20000 

20000 

20000 

20000 

15000 

1500C 

15000 

150C0 

150C0 

15000 

15000 

10000 

12000 

12000 

12000 

120C0 

12000 

8000 

8000 

8000 

8000 

80C0 


20000 

2CC00 

2CC00 

20000 

20000 

20000 

2C000 

20000 

15000 

15000 

15000 

15000 

15000 

15000 

15C00 

12000 

10000 

10000 

12000 

12000 

12000 

12000 

12000 

6000 

8000 

eooc 

8000 

8000 

8000 

4000 

4000 

5000 

5000 

5000 


UNDARIES 

EXP  PATROL 
DURATION 

WEST 

EAST 

WIN 

MAX 

0 

4000 

6.0 

7.0 

4000 

8000 

6.4 

7.4 

8000 

12000 

6.0 

7.0 

12000 

16000 

5.8 

6.8 

16000 

19000 

4.6 

7.6 

19000 

22000 

5.0 

7.0 

22000 

2800C 

5.9 

6.9 

28000 

40000 

6.5 

7.0 

0 

400C 

6.5 

7.5 

4000 

8000 

5.7 

6.7 

8000 

13000 

6.2 

7.2 

13000 

18000 

5.0 

6.5 

18000 

21000 

4.5 

7.5 

21000 

2400C 

5.3 

7.3 

24000 

3200C 

6.5 

7.5 

32000 

40000 

6.6 

6.  8 

0 

4000 

6.5 

7.5 

4000 

8000 

6.0 

7.0 

8000 

13000 

5.8 

6.  8 

13000 

18000 

5.7 

6.7 

18C00 

2200C 

4.5 

6.5 

22000 

26000 

5.0 

6.8 

26000 

32000 

6.0 

7.0 

0 

8000 

6.0 

7.0 

8000 

15000 

6.2 

7.2 

15000 

21000 

5.5 

7.0 

210C0 

24000 

4.0 

7.C 

24000 

28000 

4.5 

7.5 

28000 

4C000 

6.4 

6.6 

8000 

15000 

6.2 

7.2 

15000 

21000 

5.0 

6.5 

21000 

26C00 

4.0 

7.5 

26000 

30000 

4.0 

7.0 

30000 

40000 

6.3 

6.5 

_  TARGET  TYPE 

EXP  NS  PROBABILITIES 
OF  CPP  *  2  3  4 


8.0 
8.0 
8.0 
8.0 
10.0 
10. 0 
12.0 
12.0 
11.0 
12.0 
15.0 
18.0 
15.0 
15.0 
15.0 
10.0 
11.0 
12.0 
10.0 
11.0 
18.0 
15.0 
12.0 
8.0 
1C.0 
11.0 
22.0 
18.0 
7.0 
7.0 
8.0 
18.0 
23.0 
6.0 


0.10 

0.10 

0.10 

0.20 

0.20 

0.20 

0.20 

0.10 

0.10 

0.10 

0.10 

0.25 

0.20 

0.20 

0.10 

0.10 

0.10 

0.15 

0.10 

0.10 

0.20 

0.20 

0.20 

0.05 

0.10 

0.20 

0.20 

0.20 

0.10 

0.05 

0.20 

0.40 

0.40 

0.10 


0.20 

0.20 

0.20 

0.30 

0.40 

C.4C 

0.20 

0.15 

C.20 

0.20 

0.10 

0.25 

0.40 

C.4C 

0.15 

0.20 

0.20 

0.25 

0.10 

0.3C 

O.^o 

C.4C 

0.30 

0.35 

0.25 

0.40 

C.4C 

0.4C 

0.20 

0.35 

0.50 

0.30 

0.30 

C.20 


0.50 

0.50 

0.50 

C.4C 

0.10 

C.10 

0.50 

0.65 

0.5C 

0.50 

0.60 

C.35 

0.10 

C.10 

0.40 

0.10 

0.50 

0.30 

0.60 

0.50 

0.10 

0.10 

C.25 

0.20 

0.35 

0.20 

0.10 

C.1C 

0.10 

0.20 

0.10 

0.20 

0.20 

0.10 
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♦  REFUELING  AND  REARMING  DATA  * 


OcF-ST4TIQN  TIME  PARAMETERS  s 

MINUTES  OFF  STATION  PER  EVOLUTION  FROM  TABLE  BELOW) 


CRUI SEP  S 

DESTROYERS 

REFUEL 

REARM 

refuel 

REARM 

MIN  MAX 

MIN  MAX 

MIN  MAX 

MIN  MAX 

150  180 

240  300 

120  150 

180  24C 

INITIAL  AMMUNITION  LEVELS: 

< ROUNDS  ON  BOARD  FROM  BELOW  TABLE) 

SHIP  TYPE 

SHIP  NR  8"/55  6"/47  5H/54  5M/38 


I  8000  10000  12000  6000 


MINIMUM  ALLOWED  AMMUNITION  LEVELS  = 

8“/55  J  1000  ROS 
6«/47:  1000  ROS 

5"/54:  2000  RDS 

5H/38:  1000  RDS 


EXPECTED  AMMUNITION  EXPENDITURES: 
(ROUNOS  PER  MINUTE  PER  MISSION) 


TYPE  SHIP 


MIN  MAX 


TARGET  TYPE 

2  3  4 

MIN  MAX  MIN  MAX  MIN  MAX 


1 

2 

3 

4 


4  10 

4  11 

6  22 
4  10 


3 

3 

5 

3 


9 

10 

20 

9 


2 

2 

5 

3 


8 

8 

18 

10 


2 

2 

5 

3 


18 

10 


0)00 
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APPENDIX  F:  PRINTOUT  OF  OVERALL  STATISTICAL  RESULTS  FOR  STDF 
FEASIBILITY  SruDY  (Wi'i'H  INfUt  PARAMETERS  IAW 
APPENDIX  C) 


* *****  OVERALL  STATISTICAL  RESULTS  ****** 


TABULATED. 6ELCW  ARE  THE .MEAN  AND .THE. STANDARD.  OEVjATlON.OF  THE  MEAN 


OF  THE  INDICATED  QUANTITIES  AS  COMPUTED  FOR 
OF  THE  PROBLEM 


10  ITERATIONS 


QUANTITY 


MEAN 


STD  DEV 


TOTAL  PROBLEM  TIME  (DAYS) 

TOTAL  NUMBER  CF  MISSIONS  FIRED 


7.34 

385 


0.04 

7 


AVERAGE  WAITING  TIME  PER  MISSION  (INCL  ONLY  19.9 

MISSIONS  WHICH  WERE  REQUIRED  TO  WAIT)  (MINS) 

AVERAGE  WAITING  TIME  PER  MISSION  (ALL  MISSIONS)  5.0 

MAXIMUM  WAITING  TIME  (MINS)  95 


1.78 

0.71 

14 


AVERAGE  NUMBER  OF  MISSIONS  AWAITING  ASSIGNMENT 

MAXIMUM  NUMBER  OF  MISSIONS  AWAITING  ASSIGNMENT 

TOTAL  NUMBER  CF  MISSIONS  WHICH  WERE  REQUIRED 
TO  WAIT  FOP  ASSIGNMENT 

PERCENT  OF  MISSIONS  WHICH  WERE  REQUIRED  TO  WAIT 
FOR  ASSIGNMENT 


0.18  C. 03 

5  0.54 

95  t.71 

24.4  4.55 


AVERAGE  TIME  cROM  CALL  FOR  FIRE  (CFF)  UNTIL  FIRST 
SHOT  (INCL  ONLY  MISSIONS  REQUIRED  TO  WAIT) 

45.4 

2.12 

AVERAGE  TIME  FROM  CFF  UNTIL  FIRST  SHOT 
MISSICNS  NOT  REQUIRED  TO  WAIT) 

(INCL  ONLY 

25.7 

0.07 

AVERAGE  TIME  FROM  CF*  UNTIL  FIRST  SHOT 
(AU  MISSIONS) 

30.6 

0.78 

MAXIMUM  TIME  FROM  CFF  UNTIL  FIRST  SHOT 

132 

13.5 

AVERAGE  TRAVEL  TIME  (MINS) 

16.6 

C.  13 

AVERAGE  FIRING  TIME 

11.8 

0.27 

AVERAGE  TOTAL  SERVICE  TIME 

(TRAVEL  ♦  SETUP  ♦  FIRING  TIME) 

37.4 

0.35 
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VALUES  IN  BELCH  TABLE  ARE  MEAN  TOTAL  ROUNDS  EXPENOED  BY  EACH  SHIP 
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